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CHAPTER  1 
INTRODUCTION 

1 . 1   GENERAL 

One  of  the  most  important  characteristics  of  a  nuclear 
reactor  is  the  effect  of  temperature  on  reactivity.   As  the 
temperature  of  the  moderator,  fuel,  or  reflector  are 
individually  changed,  the  effect  on  the  reactivity  may  be 
quite  different.   Because  it  is  difficult  to  mock  up  a 
temperature  profile  in  a  low  power  reactor,  it  is  common  to 
discuss  temperature  effects  in  terms  of  the  isothermal 
temperature  coefficient  of  reactivity  where  isothermal  implies 
that  all  parts  of  the  reactor  are  at  a  constant  and  uniform 
temperature .   Comparison  of  the  measured  and  calculated  iso- 
thermal temperature  coefficient  aids  progress  toward  under- 
standing the  physical  processes  taking  place  in  the  reactor. 

This  dissertation  is  concerned  with  two  aspects  of 
thermal  effects  on  reactivity.   The  first  is  an  experimental 
study  of  the  isothermal  temperature  coefficient  of  reactivity 
of  the  Cornell  University  Zero  Power  Reactor,  a  light-water 
moderated,  2.1%  enriched,  UO^  fueled,  aluminum  clad  reactor 
having  nominal  water-to-fuel  ratios  of  1:1,  1.5:1,  2:1,  3:1, 
and  4:1.—   (See  Appendix  A  for  fuel  and  lattice  parameters.) 


—  S.  S.  Berg,  "Initial  Experiments  on  the  Cornell  University 
Zero  Power  Reactor  Cores,"  CURL-7,  p.  220  Cornell 
University  Reactor  Laboratory  (1964)  . 
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A  theoretical  prediction  is  given  in  Chapter  2,  and  the  experi- 
mental techniques  and  results  are  discussed  in  Chapter  3. 

The  second  aspect  of  this  work  singled  out  the  temperature 
dependence  of  the  disadvantage  factor  as  one  contribution  to 
the  temperature  dependent  behavior  of  the  reactor.   Chapter  4 
discusses  the  theoretical  basis  of  the  use  of  foils  to  measure 

disadvantage  factors.   Experimental  methods  and  results  are 

2/ 
discussed  in  Chapter  5  and  compared  with  predictions  by  MacVean,— 

Chapter  6  summarizes  the  results  of  both  aspects  of  this 

work  and  suggests  further  experimental  studies. 

1.2   ISOTHERMAL  TEMPERATURE  COEFFICIENT  OF  REACTIVITY 
The  isothermal  temperature  coefficient  a  is  defined  as 
dp/d0  where  p  is  excess  reactivity  and  0  is  the  temperature 
of  all  parts  of  the  core.   Isothermal  coefficients  of  reactiv- 
ity have  previously  been  measured  for  cores  similar  to  the 
Cornell  critical  assembly.   In  these  studies,  temperature 
dependent  excess  reactivities  were  determined  either  by 
measuring  periods  for  a  given  control  rod  configuration  or 
by  measuring  the  change  in  critical  position  for  a  calibrated 
control  rod. 


2/ 

-  C.  R,  MacVean,  "A  Simplified  Cell  Theory  Applied  to  the 

Calculation  of  Thermal  Neutron  Spectra  in  Light  Water 

Lattices,"  CURL-10,  p.  145  Cornell  University  Reactor 

Laboratory  (1964) . 
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The  Westinghouse  temperature  coefficient  experiments  ? 
were  made  on  aluminum  clad,  1.37o  enriched  metal  and  oxide 
fuels  of  both  0.600  and  0o387-inch  diameters.   Reactivity 
values  were  determined  by  measurement  of  periods  longer  than 
100  seconds.   Water-to-fuel  ratios  from  1.5:1  to  5:1  were 
studied  over  a  temperature  range  of  20-50°C. 

The  Yankee  cores,  also  studied  by  Westinghouse,—  had 
2.2:1  and  2.9.1  water-to-fuel  ratios  for  2„7%  enriched, 
stainless  steel  clad,  0.300-inch  diameter  UOo  fuel.   Reactiv- 
ity determination  was  by  period  measurement  with  various 
control  rod  configurations.   The  temperature  range  covered 
was  22-75 °C . 

The  Cornell  University  Zero  Power  Reactor  was  studied 
over  a  temperature  range  5-70°C  for  water-to-fuel  ratios  of 
1:1,  2:1,  and  3:1.   After  clearly  establishing  that  the  4:1 
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—  J.  R.  Brown,  et  al.,  "Reactor  Properties  of  Water  Moderated 

Slightly  Enriched  Uranium  Lattices,"  WAPD-117,  p.  117 
Westinghouse  Atomic  Power  Division  (1954)  . 

—  J.  R.  Brown,  et  al.,  "Kinetic  and  Buckling  Measurements  on 

Lattices  of  Slightly  Enriched  Uranium  or  UCU  Rods  in 
Light  Water,"  WAPD-176,  p.  39  Westinghouse  Atomic  Power 
Division  (1958)  . 

—  P.  Wo  Davison,  et_  aT.,  "Yankee  Critical  Experiments-Measure- 

ments on  Lattices  of  Stainless  Steel  Clad  Slightly  Enriched 
Uranium  Dioxide  Fuel  Rods  in  Light  Water  "  YAEC-94.  n. 
Yankee  Atomic  Electric  Company  (1959) . 
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lattice  had  a  positive  isothermal  coefficient  of  reactivity 
with  no  indication  of  going  negative,  the  study  of  this  core 
was  terminated  at  a  maximum  temperature  of  45°C.   The 
accuracy  of  the  values  of  a  for  both  the  Westinghouse  and 
Cornell  experiments  is  about  1"  0,03  4/°C,  and  this  accuracy 
is  nearly  independent  of  the  value  of  a . 

1.3   TEMPERATURE  DEPENDENT  DISADVANTAGE  FACTORS 

The  neutron  density  disadvantage  factor  5   and  the 

experimentally  measured  disadvantage  factor  5    are  defined 

by  the  ratios  5   =  N  /N  and  5    =  C  /C"  where  N  and  C   are 
J  n  exp 

respectively  the  spatially  averaged  neutron  density  and  foil 
count  rate  per  unit  volume;  superscripts  m  and  f  indicate 
respectively  the  moderator  and  fuel  regions  of  the  reactor. 
Disadvantage  factor  measurements  have  previously  been 
made  on  cores  similar  to  the  Cornell  critical  assembly .   The 
Brookhaven  measurements  y    y    y    *    were  made  with  dysprosium 


—  H.  Kouts,  "Intra-cell  Flux  Traverses,"  BNL-1783, 

Brookhaven  National  Laboratory  (1953)  . 

—  H.  Kouts,  "Intracell  Flux  Traverses  and  Thermal  Utilization 

for  1.157o  Enriched  Uranium  Rods  in  Ordinary  Water," 

BNL-1982,  Brookhaven  National  Laboratory  (1954)  . 

o  / 

—  G«  A.  Price,  "Thermal  Utilization  Measurement,"  BNL-1992, 

Brookhaven  National  Laboratory  (1954)  . 
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pellets  distributed  through  the  fuel  and  moderator  regions. 
The  pellets  were  1/16-inch  diameter,  0.007-inch  thick  foils 
made  from  a  mixture  of  Dy^Cs  with  aluminum  or  plastic. 
Numerical  integration  of  the  detailed  flux  maps  gave  average 
count  rates  for  each  region  which  were  converted  to  experi- 
mental disadvantage  factors.   An  uncertainty  of  about  2%   was 
assigned  from  data  scatter;  statistical  uncertainty  in 
counting  was  about  17o- 

The  Westinghouse  disadvantage  factor  measurements *— — 

were  made  with  uranium  sector  foils  cut  and  placed  in  the 

fuel  and  moderator  regions  of  the  reactor.   The  foil  was 

235 
highly  enriched  U    dispersed  in  aluminum  to  match  the  fuel 

235 
U    atom  density.   These  foils  were  gamma-counted  for 

fission  product  activity.   The  experimental  disadvantage 

factor  was  obtained  directly  as  the  ratio  of  count  rates  per 

unit  volume  of  the  moderator  and  fuel  sector  foils.   The 

uncertainty  in  the  measurement  was  about  3% . 


9/ 

—  H»  J.  Kouts,  £t:  al_.,  "Intracell  Flux  Traverses  and  Thermal 

Utilization,  1.0277o  Enriched  Uranium  Rods  in  Light  Water," 
BNL-1796,  Brookhaven  National  Laboratory  (1954)  . 

—  H.  C«  Honeck,  "The  Calculation  of  Thermal  Utilization  and 

Disadvantage  Factors  in  Uranium/Water  Lattices,"  Nuc . 
Sci.  Eng.  18,  49  (1964) . 
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No  previously  published  works  indicate  that  the  tempera- 
ture dependence  of  the  disadvantage  factor  has  been  studied 

experimentally.   MacVean's  recently  developed  method  of  cal- 

13/ 
culating  thermal  lattice  parameters —  gives  the  temperature 

dependence  of  the  disadvantage  factor.   Since  the  calculated 
changes  are  of  the  order  of  10   /°C  an   experimental  method 
was  required  which  was  more  accurate  than  those  previously 
used . 

Experimental  disadvantage  factors  were  measured  for 
three  water-to-fuel  ratios  over  a  temperature  range  5-65°C. 
An  integral  foil  counting  method  was  used  employing  a 
0.600  jf  0.001-inch  diameter  foil  in  the  fuel  and  an  appro- 
priately shaped  moderator  foil.   Five-mil  thick  foils  of 
dysprosium  (5.24%  dysprosium  in  aluminum),  manganese  (P-metal), 
and  vanadium  were  chosen  because  their  cross  sections  are 
nearly  1/v.   Two  counting  methods  were  used.   One  was  to 
dissolve  the  manganese  foils  in  nitric  acid  and  gamma-count 
aliquots  of  the  solutions  in  a  well-type  scintillation 
counter.   The  other  method  was  to  beta-count  the  solid  foils 
in  a  proportional  flow  counter.   Results  of  the  two  methods 
agreed  within  the  statistical  errors  which  were  about  0.7%. 


—  A.  Z.  Kranz,  "Measurements  of  Thermal  Utilization  Resonance 
Escape  Probability,  and  Fast  Fission  Factor  of  Water 
Moderated  Slightly  Enriched  Uranium  Lattices,"  WAPD-134, 
Westinghouse  Atomic  Power  Division  (1955)  . 
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— —  A.  Z,  Kranz  and  G.  G.  Smith,  "A  Second  Report  on  Measure 

ments  of  f,  p,  and  e  of  Water  Moderated  Slightly  En- 
riched Uranium  Lattices,"  WAPD-151,  Westinghouse  Atomic 
Power  Division  (1956)  . 

1 1/ 

—     Co    R.   MacVean,    footnote   2,    pp.    143-145 . 


CHAPTER  2 

THEORETICAL  CALCULATION  OF  THE  IS£EKERMAL  TEMPERATURE 
COEFFICIENT  OF  REACTIVITY 

2 . 1   INTRODUCTION 

A  calculation  of  the  isothermal  temperature  coefficient 
of  reactivity  was  made  using  two  group,  two  region  pertur- 
bation theory.   The  presentation  here  follows  closely  that 
of  Lammore  •■ — • 

Section  2.2  presents  and  discusses  the  two  group.,  two 
region  neutron  balance  equations  and  the  associated  boundary 
conditions  necessary  for  the  solution.   The  bilinearly 
weighted,  first  order  reactivity  equations  from  perturbation 
theory  are  given  for  the  core,  reflector,  and  interface 
regions  of  the  reactor.   In  these  latter  equations^  tempera 
ture  derivatives  of  the  reactor  parameters  are  weighted  by 
the  fluxes ,  adjoint  fluxes,  flux  gradients,  and  adjoint  flux 
gradients.   The  temperature  derivatives  are  discussed  in 
Sec.  2.3.   Utilization  of  these  equations  for  machine  cal- 
culations is  discussed  in  Sec.  2.4. 


J.  A.  Larrimore,  "Temperature  Coefficients  of  Reactivity 
in  Homogenized  Thermal  Nuclear  Reactors,,'1  Ph<,D,  Thesis, 
pp.  185-211,  M.I.T.  (September,  1962). 
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2„2      TWO   GROUP,    TWO  REGION    EQUATIONS 
2.2,1      Steady   State  Neutron  Balancj?_J[qua_t i^ns 

The    steady    state  neutron  balance   equations    for   the  core 
region   of   the   reactor  are: 

(DL72-Slc)    cPlc+v<ffif.s92c   =  0  (2.1) 

P2lc9lc  +    (D2cV2-22c)    cp2c   =  0  .                       (2.2) 
Analogous   equations    for   the  reflector    are: 

(Dlrv2-Slr>   *lr                         "  °  (2'3> 

^lr^lr  +   (D2r^2"22r>   *2r     "  °  (2A) 

where"  subscript  1  refers  to  the  fast  group ,  subscript  2 

refers  to  the  thermal  group,  subscript  c  denotes  the  core, 

subscript  r  denotes  the  reflector,  D  is  the  diffusion  constant 

for  a  given  group  and  region,  Z,    is  the  fast  group  removal 

cross  section,  Zr.   is  the  fission  cross  section,  S0  is  the 
'   fis  '2 

thermal  absorption  cross  section,  v  is  the  average  number 
of  neutrons  per  thermal  fission,  e    is  the  fast  fission  factor, 
^   is  the  Laplacian  differential  operator,  and  p  is  the 
resonance  escape  probability. 

The  adjoint  equations  for  the  core  and  reflector  are 
given  by 

(Dlc^2-Slc>   'Pic  +  PSlc4c   ■  °  <2-5> 

v^fis4c   +   (D2cv2--;2c>    T2c    "   °  (2-6~> 
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(DlrV2-Slr)  <p+r  +  Zlr4r  -  0  (2.7) 

(D2V2-22r)  cp+r  =  0  (2.8) 

where  the  superscript  +  indicates  the  adjoint  flux  functions 

Boundary  conditions  at  the  outer  extrapolated  boundary 
require  that  the  fast  and  slow  fluxes  and  adjoint  fluxes 
vanish .„   At  the  core-reflector  interface   the  boundary  con 

ditions  require  that  fast  and  slow  fluxes  and  adjoint  fluxes 

+     + 
match  (i0e0  cp-.   =  cp,  ,cp,   =  cp,  ,  and  similar  equations  with 

the  1  replaced  by  2),  and  that  the  fast  and  slow  currents 

and  ad -joint  currents  match  (i.e.  D.   cp.   =  D.   cp.   where  i 
J  ic   ic     ir   ir 

is  either  1  or  2,  and  the  same  expressions  are  written  for 
the  the  adjoint  fluxes), 

2.2,2   Bilinearly  Weighted  Reactivity  Relations 

Once  the  fast  and  slow  fluxes  and  adjoint  fluxes  for 
the  core  and  reflector  are  known,  they  can  be  used  in  the 
standard  two  group,  two  region,  first-order  perturbation 
theory  reactivity  equations =   Perturbations  of  the  parameters 
are  equated  to  temperature  derivatives  of  the  parameters, 
and  the  reactivity  effects  are  calculated. 

The  following  expression  is  given  by  Larrimore—'- -   for 
the  core  contribution  to  the  isothermal  temperature  co- 
efficient, a  : 
'   c 


— -  J  A.  Larrimore,  footnote  14,  p.  191 
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a      = 
c 


J    dr   cp+c    v,Sf.sCp2c\    l    Ijdr    ^c 

C  /  C 


dD-,  ,  dD0 

•  -air  V*lc  -V4c  •  t  V. 


+      ^1 
2c    '    ^lc  -3S2-  ^lc 


4 


dZ 


£•  d(pZ,    ) 

fis  ,      +         vt^    lcy 


+    <p1c  ve  -1W-    cp2c  +  cp2c  -"hq—  ^ 


+     d22c 
^c  -^0~  ^2 


3 


(2.9) 


where  dr_  is  the  volume  element,  the  subscript  c  on  the 
integral  sign  implies  integration  over  the  core  of  the 
reactor,  and  0  is  temperature. 

Similarly,  the  reflector  contribution,  a  s  is  given 
by  Larrimore—  as 

Jdr   „+.    veSfis   v2c  \     1    |  J    d£  ['V<r    • 
*  c  /  r 

dD-,      __  dD0    dZ-, 


a      = 
r 


"He" 


dZ-. 
.      +  lr 

+     9o  Jr>  CD-, 

Y2r      d9        Ylr 


'2r        ^lr  ~^0     ylr 


+     d22r 
?2r-d0"     *2r 


(2.10) 


— J.    A.    Larrimore,    footnote   14,    p.    194. 
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Finally,    Larrimore   gives    for  a.,    the   interface   effect,— 

ai  -       I  J    d£  "lc   veZfis  ^c)         •  ["Voir    -V<Plc(Dlc-Dlr) 


-Vcpt       "V?o    (D0    -D0    )    -    cpt   cp      (2,     -S-,     ) 
v  ■ 2r       v Y2c v    2c      2ry         YlrYlcv    lc      lry 

+   ^r9lc(PZlc-Slr) 

+  cp+  cp       veS  -   4rtp2c<22c-Z2r>]  .    „      ' 

-1  int 

[dV/de]  (2.11) 

where  all  the  functions  in  the  second  set  of  brackets  are 
evaluated  at  the  core-reflector  interface,  and  dV/dQ  is  the 
change  of  core  volume  with  temperature.   For  an  infinite 
cylinder,  dV/d0  is  given  by 

dV/d9  =   2irr2   B    ,/3  (2.12) 

ex  J- 

where  r  is  the  core  radius,  and  6  n  is  the  volumetric  co- 
al 

efficient  of  thermal  expansion  of  aluminum. 


— J.  A.  Larrimore,  footnote  14,  p.  209. 
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The  total  isothermal  temperature  coefficient  for  the 
reactor  is  the  sum  of  the  core,  reflector,  and  interface 
contributions : 


a  =  a  +a  -fa.  „  (2.13) 

c   r   l  v     ' 


2.3   TEMPERATURE  DERIVATIVES  OF  THE  NUCLEAR  PARAMETERS 
In  order  to  calculate  the  various  expressions  for  the 
isothermal  temperature  coefficient,  it  is  necessary  to  have 
relations  for  the  temperature  derivatives  of  the  various 
nuclear  parameters,  viz.,  v,e9p,  and  the  various  D's  and 
Z's. 

2.3.1   Reflector  Group  Constants 

The  temperature  derivative  of  the  macroscopic  cross 
section  Z  is  easily  found  by  differentiating 


Z   =  Naa  ,  (2,14) 

where  N   is  the  number  density  of  atoms  whose  microscopic 
cross  section  is  a  to  give 

(l/S)d2/d9   =    (1/NjdN    /d0  +    (l/a)da/d0    .  (2,15; 

a   a 

Since  N   is  proportional  to  the  density  d  of  the  material, 

cl 

it  follows  that 


(l/Na)dNa/d0  =  -(3  ,  (2.16) 
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where  p>  is  the  volumetric  coefficient  of  thermal  expansion 
for  the  material.   By  substitution  in  Eq.  (2.15),  the  tempera^ 
ture  derivative  of  2  becomes 

(l/Z)dZ/dQ  =  -p+(l/o)do/d0  .  (2.17) 

For  the  fast  group,  it  was  assumed  that  the  microscopic 
cross  section  does  not  change  with  temperature.   This 
assumption  leaves  only  a  density  effect  term  for  the  fast- 
group  macroscopic  cross  section  temperature  derivative.   For 
the  reflector  we  therefore  have 

(l/zlr)dZlr/d9  =  -pR  0  .  (2.18) 

The  slow  microscopic  cross  section  for  water  o^      is 

-1/2 
proportional  to  1/v;  hence  it  is  proportional  to  9     . 

Therefore, 

(l/o2r)do2r/d0  =  -0.5/0  (2. 19) 

where  0  is  in  °K.   For  the  slow  group  in  the  reflector, 
Eq.  (2.17)  therefore  becomes 

(l/22r)dS2r/d0  =  -pH  o~O.5/0  .        (2.2CT 
The  diffusion  constant  D  can  be  written 


D  =  l/(32tr)  ,  (2.21) 

where  Z        is  the  appropriate  transport  cross  section. 
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Logarithmic  differentiation  gives 

(l/D)dD/de  =  +p-dotr/de  .  (2.22) 

18/ 
Larrimore— -  gives  the  following  relationship  for  the 

logarithmic  slope  of  the  thermal  transport  cross  section 

of  water: 

(l/otr)dotr/d0  =  -n/0  ,  (2,23) 

where  n  is  given  for  20,  50,  and  100°C  as  0.473,  0.471,  and 
0,469  respectively.   Using  the  value  for  50DC,  Eq.  (2.22) 
becomes 

(l/D2r)dD2r/d0  =  PH  0+0  471/0  .         (2.24) 

Assuming  again  that  the  microscopic  cross  section  of  the 
fast  group  does  not  change  with  temperature,  Eq.  (2.22) 
for  the  fast  group  becomes 

(l/Dlr)dDlr/d9  =  6[¥r  (2.25) 

2.3.2   Core  Group  Constants 

In  the  core  region  of  the  reactor,  the  fast  group 
macroscopic  cross  section  2   can  be  written 

Zq   =    (Vf+Vmyl    (Vfif+V%)   ,  (2.26) 


18/ 

—  Jo  A,  Larrimore,  footnote  14,  pp.  83-84 
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"F  m 

where  V   and  V  are  the  volumes  of  the  fuel  and  moderator 


region  of  the  cell,  Zr   is  the  fuel  cross  section,  and 

is  the  moderator  cross  section. 

For  the  fast  group,  the  microscopic  cross  sections  are 

assumed  not  to  vary  with  temperature,  and  the  expansion  of 

the  fuel  is  taken  to  be  negligible  compared  to  the  expansion 

of  the  cladding o   The  logarithmic  temperature  derivative  of 

2,   then  becomes 
lc 

(l/2lc)dSlc/d0  =  v'V+v*)-1^/^)  |(2PA1/3)-eHJ    , 

(2.27) 
where  p.,  is  the  volumetric  coefficient  of  thermal  expansion 
of  aluminum,  the  cladding  material,   The  last  term  in 
Eq .  (2.27)  arises  because  the  density  effect  on  the  moderator 
is  the  combined  expansion  effects  of  the  water  and  aluminum 
For  the  reactor  cell,  it  is  only  the  water  cross- sectional 
area  change  per  unit  length  of  the  reactor  cell  which  con- 
tributes to  the  water  density  change*   This  area  change  is 
due  to  the  expansion  of  aluminum  structural  members.   The 
2P.-./3  term  accounts  for  this  area  change  per  unit  length. 
Using  Eq.  (2,21)  where,  for  the  fast  group, 

tr  =  (Vf+Vm)"1(VinS™r+VfZ^r)  ,  (2.28) 


Z 


it  can  be  shown  that  the  logarithmic  temperature  derivative 

of  D-,   is 
ic 
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d/Dlc)dDlc/d0  =  (Dlc/Dlr)(vm)(vf+vIT1)"1reH  0-2eA1/3) 

(2,29) 

The  temperature  derivatives  for  the  slow  group  core 
parameters  were  taken  from  calculations  by  MacVean  using  the 
computer  program  COUTH,  "a  simplified  polyenergetic  cell 
theory  (which)  is  formulated  to  determine  spatially  averaged 
energy  dependent  thermal  fluxes  in  the  moderator,  cladding, 
and  fuel  regions  within  the  unit  cell  of  a  reactor  lattice,, 
The  derived  spectra  are  then  utilized  in  the  calculations  of 

the  thermal  integral  parameters  and  aveiage  cross  sections 

19/ 
required  for  reactor  computations,"—   Calculations  at  two 

temperatures  provide  the  temperature  dependence  of  the  thermal 

parameters  used  here 

^  ° 3  . 3   Resonance  Escape  Probabili t y 

The  resonance  escape  probability  temperature  derivative 

can  be  calculated  by  the  following  scheme   Weinberg  and 

20/ 
Wigner——  give  the  following  expression  for  p(E),  the 

resonance  escape  probability  to  energy  E  from  source  energy 

E_: 


19/ 

—  C.  Ro  MacVean,  footnote  2,  p0  ii 

20/ 

— —  A.  M.  Weinberg  and  E   P  Wigner,  The  Physical  Theory  of 

Neutron  Chain  Reactors^  p..  301,   University  of  Chicago 
Press,  Chicago  (1959)  , 
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p(E)  =  exp 


L 

/ 


Za(E')[Es(E,)+Ss]"1dE,/E 


(2,30) 


where  £  is  the  average  lethargy  increase  per  collision,  2 

a 

is  the  energy  dependent  absorption  cross  section  of  the  cell, 
and  S  is  the  scattering  cross  section  of  the  cell.  Equation 
(2,31)  is  rewritten  in  the  usual  form 


-  1 


p   -    expONu(|Zs)        Rleff] 


(2,31) 


where  RI  rC   is  the  effective  resonance  integral,  and  N   is 
eff  &  u 

the  homogenized  fuel  atom  density  in  the  cell. 

Logarithmic  differentiation  with  respect  to  temperature 
results  in  the  following  expression 

(l/p)dp/de  =  -N  RI  rr'^L    )"iri/N  )dN  /dQ-(l/Z  )d   /d0  + 
r  /  f  ueft    s7       uy   u         s7   s 


(l/Rleff)dRleff 


/d9] 


(2.32") 


The  atomic  density  of  the  fuel  varies  with  temperature 
however,  since  fuel  atoms  are  lost  from  the  cell  only 
through  axial  expansion  of  the  fuel,  we  have 


(1/N  )dN  /d0  =  -6  /3 
u/      u        u 


(2,33) 


The  microscopic  scattering  cross  section  is  a  constant, 
therefore  the  scattering  cioss  section  temperature  deriva- 
tive is  given  by 
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(l/Ss)dSs/d0  =  -ps  .  (2.34) 

The  scatterer  is  water  and  is  lost  from  the  cell  by  expansion 
of  water,  but  is  gained  in  the  cell  as  the  horizontal  cross- 
sectional  area  of  the  cell  increases  due  to  aluminum 
structural  expansion .   The  density  change  in  the  moderator 
is  therefore  a  combination  of  water  and  aluminum  expansion, 
and 


(l/Ss)d2s/d0  =  -pR  0+(2/3)eM  .         (2.35) 

Using  Eqs.  (2.33)  and  (2.35),  Eq.  (2.32)  can  be  re- 
written as  follows 

(l/p)dp/d0  =  -NuRIeff (€S8)"1  [PR  0  "  T  PA1  "  7  ^u  + 


(l/RIeff)dRIeff/d9] 


(2.36) 


Table  2.1  gives  the  values  of  the  terms  in  Eq.  (2.36) 
for  the  water- to- fuel  ratios  investigated. 

2.3.4   Fast  Fission  Factor 

An  estimate  of  the  temperature  effect  on  the  fast 
fission  factor  can  be  made  using  an  empirical  relation 
for  the  fast  fission  factor  for  UOo-water  lattices  as 

U    T       ■  21/ 

given  by  Lamraore:- — 


21/ 

—  J.  A.  Larrimore,  footnote  14,  p.  145. 
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(l/Ss)dZs/d0  =  -pg  .  (2.34) 

The  scatterer  is  water  and  is  lost  from  the  cell  by  expansion 
of  water,  but  is  gained  in  the  cell  as  the  horizontal  cross- 
sectional  area  of  the  cell  increases  due  to  aluminum 
structural  expansion.   The  density  change  in  the  moderator 
is  therefore  a  combination  of  water  and  aluminum  expansion, 
and 

(l/2s)dZs/d9  =  -pR  0+(2/3)BA1  .         (2.35) 

Using  Eqs.  (2.33)  and  (2.35),  Eq .  (2.32)  can  be  re- 
written as  follows 

(l/p)dp/de  =  -NuRIeff (^8)-1  [Ph20  -  T  PA1  "  T  ^u  + 

(l/RIeff)dRIeff/d9]   . 

(2.36) 

Table  2.1  gives  the  values  of  the  terms  in  Eq .  (2.36) 
for  the  water- to- fuel  ratios  investigated. 

2.3.4   Fast  Fission  Factor 

An  estimate  of  the  temperature  effect  on  the  fast 
fission  factor  can  be  made  using  an  empirical  relation 
for  the  fast  fission  factor  for  UOo-water  lattices  as 
given  by  Lammore:-^-- 


21/ 

—  J.  A.  Larrimore,  footnote  14,  p.  145. 
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e-1  =  (0.1565) [1+0. 875(du  n  Vm/Vf)+0„288  Vs/Vf]_1    (2.37) 

where  e  is  the  fast  fission  factor,  cL  0  is  the  water  density, 

and  V   is  the  volume  of  the  structural  materials.   Logarithmic 

s   f 
differentiation,  assuming  V  /V   to  be  small,  gives 

(l/e)de/de  =  5.59(Vm/Vf)(€-l)2(€)-1[(2pH  Q/3)-0u/3]  .   (2.38) 

Using  these  values  for  the  1:1  core:  e  =  1.03,  V  /V  =  1.012, 

0  =  20°C,  pu  n  =  2.05xlO~4/°C,  and  p   =  0.3xlO~4/°C,  the 
rl^U  U 

temperature  effect  on  e  as  calculated  by  Eq.  (2.38)  is  found 

to  be  0.6x10   /°C.   For  all  other  cores  e  is  smaller  than 

2 
1.03  and  the  (e-1)   term  is  dominant.   The  fast  fission 

factor  temperature  dependence  is  therefore  found  to  be  negli- 
gible. 

2.3.5   Discussion  of  the  Temperature  Derivatives  of  the  Nuclear 

Parameters 
Temperature  derivatives  of  the  nuclear  parameters  have  two 
major  contributions:   one  which  is  directly  related  to  tempera- 
ture change  (such  as  the  variation  of  the  microscopic  cross 
sections  for  the  slow  group)  and  one  which  is  indirectly  re- 
lated to  the  temperature  change  by  the  effect  which  this  tem- 
perature change  has  on  the  thermal  expansion  of  the  materials 
in  the  reactor.   Over  the  investigated  temperature  ranges  the 
temperature  derivatives  of  the  microscopic  cross  sections  change 
by  about  10%.   Over  the  same  temperature  interval,  the  value  of 
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PA1  and  (3   are  constant,  and  dRI  ^f/d0  changes  by  about  3%„ 
Therefore,  there  are  some  terms  in  the  expression  for  the 
temperature  coefficient  of  reactivity  which  do  not  change 
appreciably  over  the  temperature  interval  studies.   Over 
this  same  temperature  range,  the  value  of  f3„  n  changes  from 
about  0,5xlO~4/°C  to  about  6,0xl0~4/°C, 

The  above  discussion  suggests  an  expression  for  a  of 
the  form 

a  =  A+BP„  ~  ,  (2  ,37) 

Calculated  values  of  A  and  B  are  tabulated  in  Table  2=3  in 
the  next  subsection.   The  experiments  described  in  Chapter 
3  showed  the  predicted  linear  dependence  although  the  cal- 
culated and  observed  values  of  A  and  B  did  not  agree  for 
all  water-to-fuel  ratios, 

2,4   NUMERICAL  CALCULATION  OF  THE  ISOTHERMAL  TEMPERATI  RE 

COEFFICIENT  OF  REACTIVITY 

22/ 
A  computer  program  was  developed  by  Osias—  which 

solves  the  neutron  balance  equations,  calculates  the  fluxes 
and  adjoint  fluxes,  and  calculates  the  isothermal  tempera- 
ture coefficient  of  reactivity  by  Eqs ,  {2  „9) -(2, 11)  for  an 


22/ 

—  Do  J.  Osias,  private  communication  (August,  1964) 
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infinite  cylinder  reactor    Calculations  were  made  at  4,  20, 
40,  60,  and  80°C.   These  calculations  represent  an  approxi- 
mation to  the  first  order  perturbation  theory  in  that  the 
flux  shapes  calculated  at  20°C  were  used  to  compute  the 
temperature  coefficient  at  all  other  temperatures.   All  para- 
meters of  the  two  group  two  region  problem  change  with 

-4  /o 
temperature  on  the  order  of  10   /  C,  so  that,  changes  in  the 

flux  shapes  with  temperature  are  expected  to  contribute  an 
error  of  less  than  \%  to  the  calculated  values  of  the  iso- 
thermal temperature  coefficient  of  reactivity, 

The  input  numbers  required  and  their  sources  are  listed 
in  Table  2.2,   Results  of  the  computer  calculation  are  given 
in  Table  2.3  for  the  five  water-to-fuel  ratios  of  the 
reactor.   The  calculated  values  of  a  are  a  linear  function 

of  PH20- 
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CHAPTER  3 

METHODS  AND  RESULTS  OF  THE  ISOTHERMAL  TEMPERATURE  COEFFICIENT 

EjgERDjENT 

3.1   GENERAL  EXPERIMENTAL  PROCEDURE 

A  typical  series  of  runs  for  the  temperature  coefficient 
experiment  was  made  by  starting  at  the  coldest  temperature  and 
progressing  to  higher  temperatures  in  steps  of  five  or  ten 
degrees.   Smaller  temperature  intervals  were  used  in  the  region 
where  the  isothermal  coefficient  changed  sign. 

Temperature  measurements  were  made  at  three  times  during 
a  run:  at  startup,  at  approximately  half  power,  and  at  maximum 
power o   The  average  of  these  temperatures  was  assigned  to  the 
run. 

All  control  rods  were  at  the  fully  withdrawn  position  and 
the  source  was  removed  from  the  core  prior  to  collecting  data 
on  the  stable  period  of  the  clean  core.   Data  were  collected 
over  approximately  three  decades  increase  in  power  level  «, 

After  the  final  temperature  measurement,  the  water  was 
removed  from  the  reactor  tank  to  the  dump  storage  tank,  heated, 
and  returned  to  the  reactor  tank  for  the  run  at  the  next  higher 
temperature . 


-27- 
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3,2   TEMPERATURE  REGULATION  AND  MEASUREMENT 
3,2,1   Water  Heating 

When  the  reactor  is  shut  down,  the  moderator  water  is  stored 
in  the  dump  storage  tank  located  just  below  the  reactor  tank ,   The 
reactor  tank  capacity  is  2,000  gallons;  that  of  the  dump  tank  is 
2 , 200  gallons.   A  steam  heater  in  the  dump  storage  tank  was  used 
to  heat  the  water  at  the  rate  of  about  one-half  degree  centigrade 
per  minute,   Since  the  system  lacks  automatic  temperature  control, 
it  was  difficult  to  obtain  or  maintain  a  given  temperature  to 
within  one  or  two  degrees „ 

32,2  Water  Chi 11 ins 


The  recirculation  pump  for  the  reactor  demineralizer  system 

circulated  the  water  from  the  dump  storage  tank  through  a 

23/ 
Worthington  water  chiller—   (model  RWW-400)  and  back  to  the 

dump  storage  tank,   The  cooling  rate  was  approximateiy  one-half 

degree  centigrade  per  hour.   Because  of  this  slow  cooling  rates 

a  series  of  experimental  runs  was  made  in  steps  of  increasing 

temperature,   The  lowest  obtainable  temperature  was  found  to  be 

about  5°C, 

3,2,3   Temperature  Measuremervt 

Four  twenty-gauge  Leeds  and  Northrup  (#58-1)  asbestos 
insulated,  fiberglass  sheathed  chromel -alumel  thermocouples 
placed  in  the  core  through  interstitial  holes  in  the  top  grid 


23/ 

—  On  loan  from  the  Department  of  Thermal  Engineering 
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plate  were  used  to  measure  the  moderator  temperature.   Sampling 
of  the  temperature  took  place  at  radii  and  depth  as  given  in 
Table  3.1. 


TABLE  3.1 

THERMOCOUPLE  POSITIONS 

Water-to- 
Fuel  Ratio 

Th 

ermocouple 
1 

Thermocouple 
2 

Th< 

2rmo couple 
3 

Th 

ermocouple 
4 

a/ 
1°1  radius- 
depth 

1.75-IV 
36 

2, 5 -IV 

18 

8=5-11 
36 

7.0- IV 
24 

2s 1  radius 
depth 

4.75-11 
30 

5.0  HI 

24 

9  .5-111 
12 

2. 25- IV 
36 

3 i 1  radius 
dep  th 

7. 5 -IV 
29 

2.  5- IV 
25 

7.5-1 
18 

5.0-III 
36 

4:1  radius 
depth 

8,75-11 

18 

£/ 

2  ,5-IV 

18 

9  .75-  III 
18 

a/ 

—  Measured  in  inches  from  the  cenCer  of  the  core.   Roman  numerals 

indicate  the  grid  plate  sector, 

b  / 

—  Measured  in  inches  down  from  the  top  of  the  upper  grid  plate 

The  horizontal  midplane  of  the  fueled  region  of  the  core  is 
at  30.5  inches. 

c/ 

—  This  thermocouple  was  located  in  the  reactor  fill  pipe  for  th] 

water-to-fuel  ratio. 


The  ice  junction  was  located  in  the  reactor  cell  next  to  the 
reactor,  and  copper  mul ticonduc tor  cable  was  used  to  bring  the 
thermocouple  voltages  into  the  control  room,   A  precision  potenti- 
ometer (Leeds  and  Northrup  model  8662)  was  used  to  measure  the 
voltages  to  "t  0.001  mV  which  corresponds  to  "t  0.02  5°C, 


Sec.  3  2  30 

3  a 2 . 4   Temperature  Uniformity 

The  time  required  for  a  step  temperature  difference  between 
the  cladding  and  the  center  of  the  fuel  rod  to  decay  to  less  than 
207o  was  computed  to  be  of  the  order  of  three  minutes.   This 
analysis  was  made  using  the  standard  Fourier  charts  for  transient 
heat  transfer  in  cylinders.   Helium  fills  the  void  between  the 
fuel  pellet  (o.da  0,600+0,003  inches)  and  the  cladding 
(id.  0.610+0  002  inches)    The  effect  of  this  void  was  included 
in  the  calculations  by  assuming  that  the  fuel  pellet  does  not 
touch  the  cladding.   This  assumption  results  in  an  overestimate 
of  the  transient  decay  time  inasmuch  as  conductive  heat  transfer 
from  the  pellet  to  the  cladding  has  been  neglected.   The  fuel 
pellet  properties  used  in  this  calculation  were: 

fuel  density  10,6  gm/cm 

thermal  conductivity     0  07  watts/cm- °C 
thermal  capacity         0  063  cal/gm-°C 

The  time  required  to  take  the  reactor  critical  is  several 
times  the  three  minute  interval  mentioned  above,   About  fifteen 
minutes  are  needed  to  pump  the  water  into  the  reactor  tank  from 
the  dump  storage  tank   The  reactor  tank  fill  pump  was  usually 
run  for  several  minutes  after  the  reactor  tank  was  full  in 
order  to  agitate  the  water  in  the  reactor  tank,   Another  eight 
minutes  are  required  to  raise  the  control  rods.   Therefore, 
normal  startup  procedure  allowed  sufficient  time  to  establish 
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thermal  equilibrium  in  the  reactor  tank,  the  maximum  observed 
temperature  difference  between  thermocouples  was  less  than  r'C» 

3  3   PERIOD  MEASUREMENT 

Following  the  establishment  of  a  good  temperature  profile, 
the  reactor  was  taken  on  a  positive  period  by  removal  of  all 
control  rods.   With  the  control  rods  withdrawn,,  the  reactor  is 
in  a  clean  configuration  because  the  fuel  followers  in  the  con- 
trol rod  elements  are  then  at  the  same  height  as  the  fueled 
section  of  the  normal  fuel  rods  = 

Reactivity  errors  resulting  from  positioning  errors 
should  be  small  because  the  differential  control  rod  worth 
is  a  minimum  at  the  fully  withdrawn  position.   The  reactivity 
effect  of  a  simultaneous  one-inch  misalignment  of  all  four 
control  rods  was  found  to  be  about  0>2(Jr„  which  is  within  the 
statistical  uncertainty  of  the  measurement.   The  specifications 
for  the  digital  indicating  system  for  control  rod  position 
call  for  position  reproducibility  to  within  0.05  inch. 
Therefore,  misalignment  of  the  control  rods  at  the  fully 
withdrawn  position  will  contribute  a  negligible  reactivity 
error,  certainly  less  than  0  *  l<j; , 

Time  dependent  flux  transients  associated  with  control 
rod  motion  or  source  removal  decay  rapidly  so  that  the  actual 

reactor  period  is  within  one  per  cent  of  the  stable  reactor 

?Zi  / 
period  after  two  periods  have  elapsed,— -   To  allow  transients 


24/ 

-1-'    Lo    J„  Templin,  "Reactor  Physics  Constants,"  ANL-5800, 

2nd  ed„,  p.  463,   Argonne  National  Laboratory  (1963) 
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to  decay,,  the  data  recording  system  was  not  started  until  the 
power  level  had  increased  beyond  that  at  criticality  by  about 
one  decade,  which  is  more  than  two  periods  in  time  , 

The  suitable  range  of  periods  was  from  25  to  100  seconds . 
Periods  less  than  25  seconds  are  too  short  to  measure  satis- 
factorily.  Lack  of  temperature,  control  of  the  water  in  the 
reactor  tank  prevented  the  use  of  periods  much  longer  than 
100  seconds  at  temperatures  much  different  from  room  tempera- 
ture because  temperature  drifts  cause  corresponding  reactivity 
drifts.   To  avoid  problems  associated  with  very  long  or  very 
short  periods,  the  excess  reactivity  of  the  clean  core  was 
kept  between  about  7(j;  and  about  25<f:  by  adding  or  removing 
peripheral  fuel  rods .   This  technique  was  chosen  rather  than 
that  of  compensating  with  control  rods  because  with  the  con- 
trol rods  fully  withdrawn,  systematic  errors  due  to  control 
rod  positioning  and  to  temperature  effects  on  control  rod 
worth  were  avoided . 

Two  BF3  pulse  counters  (Hammer- 20th  Century  31EB70G) 
provided  the  data  for  the  period  determination.   The  outputs 

of  these  counters  were  recorded  automatically  by  a  Hamner 

25/ 
Spectrometer,  a  combination  scaler-timer-printer .■ —   The 

data  recording  cycle  was  adjusted  so  that  several  data  points 

were  recorded  during  the  time  interval  of  one  period.   The 


— '  S.  S.  Berg,  footnote  1,  p.  104 
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detectors  themselves  were  positioned  in  the  reactor  instrument 
thimbles  so  that  saturation  did  not  occur  for  power  levels  used. 

3.4   DATA  REDUCTION 

Period  measurement  raw  data  were  analyzed  by  the  computer 

26  / 
code  RHOFIT,— -  a  least  squares  fit  to  an  exponential  function. 

The  computed  period  is  then  fed  into  a  subroutine  which  cal- 
culates the  excess  reactivity  by  the  inhour  equation 

6 
Pex^eff  -   A^eff^"1      2    (Pi/Pef  f  >  <1+AiT) "'  (3^> 


i=l 


where 


A  is  the  neutron  generation  time 

B.  is  the  importance  weighted  effective  delayed  neutron 
yield  for  the  i    group 

t"V» 

A.  is  the  decay  constant  for  the  precursor  of  the  i 

group 

B  rr   is  the  effective  delayed  neutron  fraction 

T  is  the  stable  period  of  the  reactor 

p   /B  .cc  is  the  excess  reactivity  in  dollars. 
r ex  err  J 

Table  3.2  lists  the  values  of  the  parameters  used  in  the  inhour 
equation  as  well  as  selected  values  of  p/8  rr   and  T  for  the 
four  water-to-fuel  ratios  used. 


9  (\  I 

—  W.  E.  Schilling  and  D.  H.  Bryce,  private  communication 
(March,  1964) . 
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TABLE  3.2 


REACTIVITY  PARAMETERS- 


W/F 


1:1 


2:1 


3:1 


4:1 


?!<%) 

0.0397 

0.0317 

0.0310 

0,0321 

P2(7o) 

0.2206 

0.2125 

0.2142 

0,2147 

P3(%) 

0.2023 

0.1933 

0.1942 

01946 

P4(%) 

0.3553 

0.3955 

0.3945 

0,3945 

P5(%) 

0.1342 

0.1229 

0.1216 

01207 

P6(%) 

0.0480 

0.0440 

0.0437 

0,0435 

Peff^%\ 

0.7943 

0.8630 

0.8330 

0,7955 

A -j  (sec   ) 

0.0124 

0.0124 

0*0124 

0.0124 

A^^ec"  ) 

0.0305 

0.0305 

0.0305 

0.0305 

Ao(sec   ) 

0.1111 

0.1111 

0.1111 

0.1111 

\    t        "1\ 

A,(sec   ) 

0.301 

0.301 

0.301 

0,301 

■v  /    -1\ 
Ac-(sec   ) 

1.14 

1.14 

1.14 

1.14 

A,- (sec   ) 

3.01 

3.01 

3.01 

3.01 

A/B  ..(sec) 

0.004041 

0.005411 

0.007155 

0.008976 

Reactivity  (in  $)  for  period  T  (in  seconds) 


T(sec) 


20 

0.2889 

0.2799 

0.2812 

0.2819 

30 

0.2307 

0.2223 

0.2234 

0.2240 

40 

0.1937 

0.1859 

0.1868 

0.1874 

50 

0.1676 

0.1604 

0.1612 

0.1617 

60 

0.1480 

0.1414 

0.1421 

0.1425 

70 

0.1327 

0.1265 

0.1272 

0.1276 

80 

0.1204 

0.1146 

0.1152 

0.1156 

90 

0.1102 

0.1048 

0.1054 

0.1057 

100 

0.1017 

0.0967 

0.0971 

0.0974 

-  S.  S.  Berg,  footnote  1,  pp.  157,  210,  232  (corrected  by 
D.  H.  Bryce) . 
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The  difference  in  reactivity  for  two  runs  was  divided  by  the 
difference  in  temperature  to  obtain  values  of  the  isothermal 
coefficient  a.   The  average  temperature  0  and  the  corresponding 
value  of  the  volumetric  coefficient  of  expansion  of  water  6„  ^ 
for  the  two  runs  were  assigned  to  that  value  of  a.  L  ,  was 

determined  by  differencing  a  table  of  temperature  dependent 

27/ 
volumes  of  water.—   Figure  B.l  is  a  plot  of  &„  n  vs  0.   Tables 

H2° 

B.l  -  B,4  in  Appendix  B  give  the  temperature,  period,  excess 

reactivity,  and  number  of  fuel  rods  for  each  run  and  for  all 

water- to- fuel  ratios   Values  of  a  and  FH  0  given  in  the  same 

tables  were  used  in  another  computer  code  to  least  squares  fit 

the  values  of  a  to  a  linear  function  of  Fu  ^ > 

ri  ^U 

3.5   EXPERIMENTAL  RESULTS 
3.5.1   Temperature  Dependent  Excess  Reac tivi_ty_ 

Figures  3.1  through  3.4  show  the  temperature  dependent 
excess  reactivities  of  the  cores  investigated.   Uncertainties 
in  the  values  of  p  are  within  the  size  of  the  points  of  the 
graphs .   Nominal  water-to-fuel  ratios  are  designated  X:l-2 
(read  X  to  one  dash  two)  where  X  is  the  water  relative  volume 
and  the  -2  indicates  the  core  configuration.   All  the  -2  cores 
were  set  up  in  the  reactor  with  the  removable  unit  cell  at  the 


27/ 

— —  C=  D.  Hodgman  (Editor-in-Chief),  Handbook  of  Chemistry  and 

Physics,  42nd  ed ,  ,  p.  2143,   Chemical  Rubber  Publishing 

Co.,  Cleveland,  Ohio  (1960), 
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center  of  the  lattice.   As  a  result,  all  -2  cores  are  off  center 
with  respect  to  the  reactor  tank.   A  minimum  distance  of  more 
than  two  feet  separated  the  closest  fuel  rods  from  the  grid  plate 
support  structure.   The  temperature  dependent  behavior  of  the 
reactor  is  not  expected  to  be  affected  by  the  eccentricity  of 
the  core  in  the  reactor  tank. 

3.5.2   Isothermal  Temperature  Coefficient  of  Reactivity 

The  isothermal  reactivity  coefficient  for  each  of  the  four 
water-to-fuel  ratios  is  shown  in  Figure  3.6  as  a  function  of 
temperature  and  in  Figure  3.7  as  a  function  of  the  volumetric 
coefficient  of  expansion  of  water.   Figures  3.8-3.11  give  these 
values  of  a  vs  ]3„  n  in  greater  detail  for  each  core.   The 
straight  lines  in  Figures  3.7-3.11  represent  the  equation 

a  =  A+BFH  0  (3.2) 

where  A  and  B  were  determined  by  the  method  of  least  squares. 
Fitted  values  of  A  and  B  are  given  in  Table  3.3. 

3.6   DISCUSSION  OF  RESULTS 
3.6.1   Dissolved  Air  in  the  Water 

Some  previous  temperature  coefficient  work  was  plagued 

28  2  9/ 
with  the  formation  of  air  bubbles  on  the  fuel  rods."—* — - 


28  / 

—  J.  R,  Brown,  et  al.,  footnote  4,  p.  40. 

29/ 

—  P.  W.  Davison,  et  al . ,  footnote  5,  p.  69 
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Figure  3.4:   Excess  Reactivity  vs  Temperature  for  Various 
Fuel  Rod  Loadings,  N;  4:1-2. 
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vs  Volumetric  Coefficient  of  Expansion  of  Vater;  1:1 
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vs  Volumetric  Coefficient  of  Expansion  of  Water;  4:1-2. 
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TABLE  3o3 
FITTED  LEAST  SQUARES  VALUES  OF  A  AND  B  FOR  THE  EQUATION 

A+B\o 


W/F 


A(in  4/°C) 


B(in  10" 


Temperature 
Interval 
(°C) 


1:1 

2:1 
3:1 
4:1 


0.29+0.02 
0.53+0,02 
0.82+0,02 
1.24+0  03 


0.435+0,005 

■0,367+0.006 

0,268+0.005 

■0.138+0.011 


5  -  70 

5  -  78 

5  -  62 

7  -  45 


These  small  air  voids  were  thought  to  be  responsible  for  some 
unexpected  results  in  the  measurement  of  the  temperature 
coefficient.   Two  precautions  were  used  by  previous  investi- 
gators  One  was  to  use  deaerated  water  to  prevent  the  problem 
from  arising.   The  other  solution  was  to  start  at  maximum 
temperature  and  cool  the  reactor.   Air  solubility  being  a 
maximum  at  low  temperatures,  by  cooling  the  water  from  high 
temperatures  no  dissolved  air  would  be  released  , 

Two  precautions  were  taken  in  the  present  experiment  to 
prevent  the  formation  of  air  bubbles  on  the  fuel  rods-   First, 
the  fuel  rods  were  cleaned  with  acetone  to  remove  grease  and 
fingerprints  which  might  tend  to  collect  air  bubbles-   Second, 
water  pumping  was  continued  after  the  water  level  reached  the 
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mouth  of  the  overflow  pipe  leading  back  to  the  dump  storage  tank 
The  resulting  circulation  of  water  in  the  reactor  tank  tended  to 
flush  away  any  air  bubbles  (and  also  speed  up  the  approach  to  the 
equilibrium  temperature) . 

Systematic  errors  due  to  variation  in  dissolved  air  content 
of  the  water  were  estimated  to  be  negligible,   In  the  Cornell 
Zero  Power  Reactor,  all  water  heating  and  cooling  takes  place 
external  to  the  reactor  tank.   The  only  water  temperature  change 
occurring  in  the  reactor  tank  is  a  small  one  (1°C  or  less)  which 
is  brought  about  by  heating  or  cooling  of  the  fuel  rods  and  the 
reactor  structure  by  the  water  as  it  is  pumped  into  the  reactor 
tank, 

An  agitator  is  installed  in  the  dump  storage  tank  to  insure 
uniform  properties  of  the  water  when  pumped  to  the  reactor  tank 
This  mixing  action  aerates  the  water;  therefore  it  was  assumed 
that  the  water  at  each  temperature  was  a  saturated  solution  of 
air  and  water . 

"Between  the  temperatures  of  5°C  and  8°C,  the  density  of 

water  saturated  with  air  was  found  to  be  0,000003^  gm/ml  less 

30/ 
than  the  density  of  air-free  water," —   At  temperatures  above 

this  range,  the  solubility  of  air  in  water  is  less  and  therefore 


30/ 

—  E.  W,  Washburn  (Editor-in-Chief),  International  Critical 

Tables  of  Numerical  Data,  Physics,  Chemistry  and  Technology, 

1st  ed„,  p.  26.   McGraw-Hill  Book  Company,  Inc „ ,  New  York 

(1928)  . 
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the  density  effect  is  presumably  even  smaller.   Density  changes 
due  to  thermal  expansion  are  large  by  comparison  (of  the  order 
10"4/°C)o 

Under  the  assumption  that  dissolved  air  is  added  or  removed 
in  the  dump  storage  tank,  the  change  in  solubility  of  air  in  the 
water  results  in  a  density  change  which  is  of  no  consequence  for 
th i  s  e xp  er imen  t . 

3  . 6 . 2   Temperature  Dependent  Excess  Reactivity 

The  excess  reactivity  has  been  determined  for  four  water- 
to-fuel  ratios  over  the  temperature  intervals  given  in  Table  3  3 
Figures  3.1-3.4  give  temperature-dependent  excess  reactivities 
for  various  fuel  rod  loadings  with  sufficient  accuracy  to  plan 
further  temperature  dependence  studies.  Four  thermocouples 
with  a  total  worth  of  about  -2cf:  were  in  the  lattices  for  these 
measurements  . 

For  all  water-to-fuel  ratios,  the  excess  reactivity  curves 
have  an  initial  positive  slope.   The  temperature  at  which  the 
maximum  excess  reactivity  was  observed  was  8.8,  15.5,  and  30.2 
for  the  1:1,  2:1,  and  3il   water-to-fuel  ratios  respectively. 
For  the  4:1  core  no  maximum  was  reached;  a  curve  of  excess 
reactivity  vs  temperature  for  this  core  is  shown  in  Figure 
3oll»   Experimental  data  are  shown  up  to  about  45°C;  for  higher 
temperatures  the  plotted  points  are  from  a  linear  extrapolation 
of  a  vs  (3.   For  an  increase  in  isothermal  temperature  from  20 'C 
to  100°C,  the  reactivity  added  would  be  354  >   Thus,  for  normal 


Sec,  3.6  51 

room  temperature  operating  conditions,  a  temperature  increase  to 

100°C  will  not,  in  itself,  lead  to  prompt  criticality ,   Since 

31/ 

the  Safeguards  Analysis—  was  based  on  the  4:1  core  and  assumed 

prompt  criticality,  there  is  no  significant  change  in  the 
analysis  brought  about  by  the  determination  of  a  positive  iso- 
thermal temperature  coefficient  of  reactivity  for  this  core 

3  . 6  o  3   Isothermal  Temperature  Coefficient  of  Reactivity 

The  isothermal  temperature  coefficient  of  reactivity  has 
been  measured  for  water-to-fuel  ratios  of  1:1,  2:1,  3:1,  and 
4:1  over  temperature  intervals  listed  in  Table  3,3.   The 
accuracy  of  the  individual  measurements  of  a  is  about  0.03- 
0 ,06(j:/oC  „   The  accuracy  determined  for  a  from  the  least  squares 
fitted  curve  is  0.02(|:/oC  for  three  lattices,  and  0,03^/°C  for 
the  4:1  lattice. 

Within  the  accuracy  of  the  measurements,  A  and  B  are 
different  constants  for  each  water-to-fuel  ratio.   Theoretical 
and  experimental  values  of  A  and  -B  are  plotted  vs  water-to- 
fuel  ratio  in  Figure  3.12.   The  experimental  curves  shown 
merely  connect  the  measured  points  smoothly  and  have  no 
theoretical  basis.   However,  they  permit  estimates  by  inter- 
polation for  the  1,5:1  water- to- fuel  ratio:  A  =  0  .40  JO, 02  in 
4/°C  and  B  =  -0,405+0.006  in  104^ 


31  / 

—'    S,  S,  Berg,  "Final  Safeguards  Report  to  the  U,  S   Atomic 

Energy  Commission  for  the  Cornell  University  Zero  Power 

Reactor,"  CURL-4,  pu  7-3,   Cornell  University  Reactor 

Laboratory  (1962) . 
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Theory  predicts  and  experiment  verifies  a  linear  relation 
ship  between  a  and  ]3„  n;  however,  Figure  3,12  shows  disagree- 
ment between  the  measured  and  calculated  values  of  A  and  B. 
This  is  a  well-defined  experiment.   All  parts  of  the  reactor 
were  at  the  same  temperature.   The  control  rods  were  fully 
withdrawn,  and  the  fuel  followers  provided  a  clean  core. 
For  each  temperature,  a  stable  reactor  period  was  measured, 
There  were  no  large  perturbations  introduced  by  experimental 
devices.   The  experiment  does  not  interfere  with  itself,  as 
do  flux  measurements  using  foils.   The  conclusion  therefore 
must  be  that  the  theory  fails  to  predict  the  proper  results, 
and  further  refinements  in  the  calculations  are  necessary 

The  only  available  previously  published  data  which  were 
complete  enough  and  covered  a  large  enough  temperature  range 
to  permit  useful  plots  of  the  isothermal  temperature 
coefficient  of  reactivity  as  a  function  of  temperature,  and 
as  a  function  of  the  volumetric  coefficient  of  expansion  of 

water,  were  for  the  SPERT-I  reactor,  a  highly  enriched  pool 

32/ 
reactor,—  The  temperature  range  covered  was  15-97°C. 

Figure  3,13  shows  both  a  vs  0,  and  a  vs  Fu  r\   f°r  the 
SPERT-I  reactor,   The  temperature  coefficient  of  reactivity 
is  seen  to  be  a  linear  function  of  the  volumetric  co- 
efficient of  expansion  of  water  for  this  reactor  also 
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3.7   CONCLUSIONS 

The  excess  reactivity  has  been  measured  as  a  function 
of  temperature  for  water-to-fuel  ratios  of  1:1,  2:1,  3:1, 
and  4:1  with  sufficient  accuracy  to  allow  planning  of  further 
temperature  dependence  experiments . 

The  isothermal  temperature  coefficient  of  reactivity 
has  been  measured  for  the  same  water-to-fuel  ratios  and  was 
found  to  be  a  linear  function  of  the  volumetric  coefficient 
of  expansion  of  water  for  all  lattices.   Furthermore,  the 
results  of  SPERT-I  also  exhibit  this  linearity  of  a  vs  |3   ~ . 
Theory  disagrees  with  measured  values  of  A  and  B,  indicating 
a  need  for  more  sophisticated  calculations  for  the  isothermal 
temperature  coefficient  of  reactivity. 


CHAPTER  4 
DISADVANTAGE  FACTOR  MEASUREMENT  THEORY 
4.1  THEORETICAL  DISADVANTAGE  FACTORS 

Theoretical  thermal  disadvantage  factors  are  calculated 
by  taking  ratios  of  spatially  averaged  thermal  neutron  den- 
sities, which  include  neutrons  of  all  energies  below  the 
thermal  neutron  cutoff  energy  E   (or  below  the  corresponding 

velocity  v  ) . 

33/ 
Honeck —  gives  the  following  definitions  for  the  thermal 

neutron  density,  N(r_)  ;  thermal  neutron  density  disadvantage 

factor,  5  ;  and  the  thermal  flux  disadvantage  factor,  5  : 

?      n  cp 

v 

r 

N(r)  =  /    N(r,v)dv  (4.1) 

0 

(1/Vm)  J   dr  Nm(r) 
5   = 22^ , .  (4.2) 


n 


(1/Vf)    f     dr  Nf (r) 


where 


fuel 


5   =  N111/^  (4.3) 

n 

5^  =  ^7?f  (4.4) 


V  V 

v=  /     vN(v)dv/  J    N(v)dv  (4.6) 

0  0 


— /  H.  C.  Honeck,  footnote  10,  p.  50 

-56- 
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in  the  fuel  and  moderator  regions,  and  where: 
r_  is  the  position  vector; 

-F 

V  ,  V   are  the  moderator  and  fuel  volumes; 

-f- 
N  ,  N   are  the  spatially  averaged  thermal  neutron  densities 

in  the  moderator  and  fuel  regions; 

rn    — f 

cp  ,  cp   are  the  spatially  averaged  total  thermal  fluxes  in 

the  moderator  and  fuel  regions;  and 
v  ,  v   are  the  velocities  averaged  according  to  Eq.  (4o6) 

in  the  moderator  and  fuel  regions . 

34/ 
MacVean— —  used  similar  equations  to  calculate  the 

thermal  neutron  densities,  fluxes,  average  velocities, 

density  disadvantage  factors  and  flux  disadvantage  factors 

for  the  five  water-to-fuel  ratios  of  the  Cornell  University 

Zero  Power  Reactor.   His  results  will  be  given  later  in 

this  chapter  as  their  use  is  required. 

4.2   EXPERIMENTAL  DISADVANTAGE  FACTORS 
The  experimentally  measured  quantity  6    which  is  re- 
lated to  the  neutron  density  disadvantage  factor  is  here 
called  the  dysprosium,  manganese,  or  vanadium  disadvantage 
factor  depending  on  the  foils  used.   This  quantity  is  defined 
by 


34/ 

—  C.  R.  MacVean,  footnote  2,  pp.  143-145. 
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count  rate  per  unit  volume  in  the  moderator  foil 

5     =  , _ , , — __ — _ _ ___  . 

P     count  rate  per  unit  volume  in  the  fuel  foil 

(4.7) 

For  a  foil,  the  count  rate  per  unit  volume  C  is  given  by 

oo 


hi 


C  =  e  a  s  (1/V)  K(t)  \  d r  f   dE  5a(E)cp(E,r) 


0 

(4.8) 


where 


e  is  the  counter  efficiency; 
a  is  the  foil  self-absorption  correction; 
s  is  the  foil  self-shielding  correction; 
V  is  the  foil  volume  (integration  over  r_  is  over  the 
foil  volume) ; 
2  (E)  is  the  energy  dependent  activation  cross  section; 
cp(E,r_)  is  the  energy  and  position  dependent  flux;  and 

K(t)  is  a  time  dependent  function  which  takes  into  account 
the  irradiation  time,  the  time  since  reactor  shut- 
down, and  the  counting  time.   For  any  given  dis- 
advantage factor  measurement,  K(t)  will  be  the  same 
for  both  moderator  and  fuel  foils  and  will  cancel 
out  of  the  expression  for  the  disadvantage  factor. 
In  practice  it  need  not  be  calculated. 
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The  counter  efficiency,  e,  will  be  different  for  the 
differently  shaped  foils.   This  effect  and  an  experimentally 
determined  correction  is  discussed  in  the  next  chapter.   For 
purposes  of  this  chapter,  e  is  assumed  to  cancel. 

For  small  source-to-counter  distances,  the  self- 
absorption  factor,  a,  is  a  function  of  the  source  thickness 

35/ 
only. —   Any  sensitivity  to  geometric  shape  is  determined 

experimentally  and  included  in  the  counter  efficiency 

correction  factor  discussed  in  the  next  chapter.   For  purposes 

of  this  chapter,  then,  the  self-absorption  factors  will  cancel 

in  the  expression  for  the  disadvantage  factor. 

The  self-shielding  factor,  s,  may  vary  slightly  from  the 
fuel  to  moderator  foils,  but  this  variation  is  neglected  for 
this  work  and  the  self-shielding  factors  are  assumed  to  can- 
cel in  the  expression  for  the  disadvantage  factor.   Experi- 
mental determination  of  the  self-shielding  effect  could  be 
accomplished  by  measuring  5    as  a  function  of  foil  thickness 
and  extrapolating  to  zero  thickness. 

If  the  energy  range  is  broken  up  into  the  thermal  region 
(E  varying  from  0  to  E  ),  and  the  fast  region  (E  varying  from 
E   to  cd  ) ,  the  following  expression  can  be  written  for  a  cal- 


culated value  5   -.   of  the  experimental  quantity 

CalC 


b 
exp 


35/ 

—  W.  J.  Price,  Nuclear  Radiation  Detection,  p.  132.   McGraw- 
Hill  Book  Company,  Inc.,  New  York,  Toronto,  and  London 
(1958). 
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E  CD 

c 
e   a   s  Ut)^)'1  Jdrj^f     dE2a(E)cpm(E,  r)   +  /  dE-a(E)cpm(E,r)j 

° fc 

0         i  =    — p ■ ■ — 

calc  E  co 

c 

e   a   s   K(t)(Vf)"1Jdr  I    J       dESa(E)cpf  (E,  r)    +  J    dE^a(E)cpf  (E,r^  I 

0  Ec 

(4.9) 

On  canceling  e,  a,  s,  K(t),  and  on  using  the  following  relation 

for  the  energy  dependent,  spatially  averaged  flux, 


cp(E)  =  V"1  Jdrcp(E,r)  ,  (4.10) 


Equation  (4.9)  can  be  written  as 


E„  co 

.c 


/   dESa(E)^n(E)  +  /  dESa(E)9m(E) 

0  E 

5   ,   =  = ^ .      (4,11) 

calc      E  co  v     ; 

c 

/   dEZa(E)?f(E)  +  /   dESa(E)9m(E) 

0  E 

c 


4.2.1   1/v  Thermal  Cross  Section  Foils 

For  a  foil  with  a  1/v  thermal  activation  cross  section, 


2  (E)  can  be  written 
a v 


2  (E)  =  Z   v  /v  (4.12) 

ax      o  o  ' 


where  the  subscript  o  indicates  a  reference  energy,  and  2  v 
is  a  constant.   The  calculated  disadvantage  factor  then 
becomes 
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oo 


f      dESovQv- VCE)    +  f       dEZa(E)^pm(E) 


0 


calc,l/v  E 


oo 


/      dE2^v^v"1cpf(E)    +    /      dE2„(E)cpf(E 


o   o 


/ 


0 


(4.13) 


Using   the   relations   cp/v  =  N,    cp  =  C./E    for  E   >  E      (where  C. 

J  c  j 

is    a  constant),    and 


v 


c  c 

J       N(v)dv   =   /      N(E)dE    , 

0  0 


(4.14) 


it    follows    that 


CD 


S   v  N™  +  Cm 
o   o 


™     /       dE   E"1Sa(E) 


calc , 1/v 


oo 


2   v  Nf  +  C?  dEE    LS    (E) 

o   o  j     *  av 


(4.15) 


By  simple  manipulations  it  can  be  shown  that 

,m     „  -r  /  / 


5      i       -\  i      =6 
calc, 1/v  n 


1  +  C1"  Rl/(voaoN111) 


± 


1  +  C^   Rl/(voooNf)_I 


=    5   W 
n 


(4,16) 


Sec.  4.2  62 

where 

5   is  defined  by  Eq.  (4.3); 

m      -F 

C-,  C  are  the  epi thermal  flux  proportionality  constants; 

W  is  defined  by  the  expression  in  the  brackets; 

R.I  is  the  infinite-dilution  resonance  integral  for 

the  foil; 

v   is  the  thermal  reference  velocity  (2200  m/sec);  and 

o   is  the  activation  cross  section  of  the  foil  at  the 
o 

energy  corresponding  to  v  . 

If  the  second  term  in  both  the  denominator  and  numerator 
of  Eq.  (4ol6)  were  small  compared  to  1,  5   -,   -,  /   would  be 
equal  to  5  . 

The  resonance  integral  for  the  foil  is  defined  by 

oo 


/ 


E_1aa(E)dE  .  (4.17) 


E 


For  the  experimental  values  of  RI  used  in  the  present  cal- 
culations, E  was  0,4  eV.   However,  the  values  of  N,  6  , 
'   c  '  '   n' 

and  6   calculated  by  MacVean  are  for  a  cutoff  energy  of 

1.6  eV.   Corrections  for  the  different  values  of  E  were 

c 

made  according  to  the  following  scheme,, 
Let  N-,  /-  be  given  by 

1.6 
N1<6  =  J    dE  N(E)  .  (4.18) 

0 
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A  similar  expression  defines  Nq  ,  .   Then 

1.6 
NQ  4  =  N1#6  -  J       dE  N(E)  .  (4.19) 

0.4 

In  the  region  0.4  :  E    1.6,  let 

N(E)  =  cp(E)/v(E) 

=  C.E'^CE)"1  (4.20) 

where  C.  is  selected  to  be  that  value  given  by  Ecp(E)  from  the 
COUTH  spectra—7  at  0.55  eV.  For  v(E)  in  m/sec,  and  E  in  eV, 
the  relationship  between  v  and  E  can  be  written 

0.723xl0"6  v(E)  =  E1/2  .  (4.21) 

Using  Eqs .  (4.20)  and  (4.21),  Nq  /  is  seen  to  be 

^0.4  =  N1.6-1.^2xl0-6  C.  (4.22) 

where  N  is  in  n/cm  and  C.  is  in  n/cm  -sec. 

Table  4.1  gives  the  values  of  o  ,  v  ,  and  RI  for  the 
to  o'   o* 

three  foil  materials  used. 

Values  of  W,  C.,  Nq  ,,  and  the  calculated  &   -,   -,  / 
are  listed  in  Table  4.2  for  the  foils  used. 

4o2o2   Foils  with  a  Non-l/v  Thermal  Cross  Section 

For  a  foil  whose  thermal  cross  section  deviates  from  1/v, 

5   ,   can  be  computed  from  Eq.  (4.11)  by  integrating  the 
ca  jlc 

energy  dependent  cross  sections  over  the  spectra  of  the  fuel 
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TABLE  4.1 


164 


RESONANCE  INTEGRALS  AND  2200  M/SEC  CROSS  SECTIONS  FOR  Dy 

Mn,  and  V 


Dy164       Mn        V 

RI  (in  barns)  482^     14. 3-/      2-/ 

aQ  (in  barns)  260C£/     13. 4^    4.5^ 


a/ 

—  J.  L.  Crandall  in  "Reactor  Physics  Constants,"  2nd  ed. 

(L.  J.  Templin,  ed . )  ANL-5800,  p.  164,  Argonne 
National  Laboratory  (1963) . 

h  / 

—  D.  J.  Hughes,  Pile  Neutron  Research,  p.  139.   Addison- 

Wesley  Publishing  Co.,  Inc.,  Cambridge,  Mass.,  1953. 

c/ 

—  D.  J.  Hughes  and  J.  A.  Harvey,  "Neutron  Cross  Sections," 

BNL-325,  pp.  7-18,  Brookhaven  National  Laboratory 
(1955). 
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TABLE  4.2 

THERMAL  PARAMETERS  FOR  CALCULATING  DISADVANTAGE  FACTORS  AT  20 °C 

FOR  1/v  THERMAL  CROSS  SECTIONS 


W/F 

1:1 

2:1 

3:1 

mod    fuel 

mod    fuel 

mod    fuel 

a/      2 
C.  —   (n/cm  -sec) 

.7993  .7708 

.8274  .7955 

.8371  .8028 

NL  6  -^    (in  10"6n/cm3)   14.58  11.65   28.62  21.18   41.50  29.19 
NQ  4  -1    (in  10"6n/cm3)   13.67  10.77    27.68  20.27   40.54  28.27 

5n,0.4 
Dy 

6Dy,l/v 
Mn 
°Mn,l/v 
WV 
6V,l/v 

-  C  =  Ecp(E)  @    0.55  eV  from  COUTH  Spectra:  C.  R.  MacVean, 

footnote  2,  pp.  136-139. 

-  W.  E.  Schilling,  private  communication  (August,  1964). 

-  Calculated  from  Eq.  (4.22). 


1.269 

1.366 

1.434 

0.990 

0.992 

0.993 

1.256 

1.355 

1.424 

0.953 

0.962 

0.967 

1.209 

1.314 

1.387 

0.976 

0.982 

0.986 

1.238 

1.341 

1.414 
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and  moderator.   The  expression  for  5   -,   can  be  simplified 

CalC 

somewhat  by  defining  appropriate  average  cross  sections  and 
utilizing  the  resonance  integral  information  above. 

The  spatially  averaged  total  thermal  neutron  flux  for 


both  the  moderator  and  fuel  regions  is  given  by 

Ec 
9  =  J       9(E)dE  . 

0 

Average  thermal  cross  sections  are  defined  by 


(4.17) 


=  /   dEaa(E)cp(E)/j"   dEcp(E) 


0 


(4.18) 


where  o"  for  the  moderator  and  fuel  regions  are  integrated 

over  the  respective  spectra. 

It  can  be  shown  that  application  of  the  above  expressions 

gives  for  5   i 
b  calc 


m 


calc 


=    6 


a 
a 


,m 


m  -m. 


1  +  Cm  Rl/(al"  ?") 
J . 2 


1  +  Cf5    Rl/(?  9f) 


m/"T> 


VvO  x 


(4.19) 


where  X  is  defined  by  the  expression  in  brackets.   For  a  1/v 
thermal  absorber, 
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a  = 
a 


E  E 

/c  ~c 

dE  aovocp(E)v-1/ J   dEcp(E) 
0  0 


=  a  v  /v"  (4.20) 


o  o 


and  Eq.  (4.19)  reduces  to  Eq.  (4.16)  as  it  should. 

Both  manganese  and  vanadium  thermal  activation  cross 
sections  are  very  nearly  1/v.   The  thermal  activation  cross 

section  for  dysprosium,  however,  does  vary  from  1/v.   The 

37/ 
dysprosium  cross  sections  were  taken  from  Sher,  et  al . — 

These  dysprosium  cross  sections  and  1/v  cross  sections  for 

manganese  and  vanadium  were  averaged  over  the  COUTH  spectra 

for  water-to-fuel  ratios  of  1:1,  2:1,  and  3:1.   The  values 

of  these  averaged  cross  sections,  values  of  X,  qL  /  (total 

thermal  flux  for  a  cutoff  energy  0.4  eV) ,  and  the  various 

disadvantage  factors  are  given  in  Table  4.3. 

4.3   TEMPERATURE  DEPENDENCE  OF  THE  DISADVANTAGE  FACTOR 
MacVean  calculated  the  thermal  neutron  density  and 
flux  disadvantage  factors  according  to  the  equations  of 
Section  4.1.   By  calculating  these  thermal  properties  of 


—  R.  Sher,  et  al . ,  "Low  Energy  Cross  Sections  of  Dy    ," 
Nuc.  Sci.  Eng.  11,  369-376  (1961). 
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TABLE  4.3 

THERMAL  PARAMETERS  FOR  CALCULATING  DISADVANTAGE  FACTORS  AT  20°C 
FOR  NON-l/v  THERMAL  CROSS  SECTIONS 


W/F 

1:1 

2:1 

3:1 

mod    fuel 

mod    fuel 

mod     fuel 

H1.6  -      (n/cm3) 

14.58   11.65 

28.62   21.18 

41.50   29.19 

Vj_  6  -/  (m/sec)  3667  3971  3134  3446    2945   3265 

9-l  6  (n/cm2-sec)  5.3465  4.6262  8.9695  7.2986  12.2218  9.5305 

C.  -     (n/cm2-sec)  0.7993  0.7708  0.8274  0.7955  0.8371  0.8028 

~     ~     -     (barns)  1735  1610  1937  1765    2019    1820 


a,Dy 
a,Mn 
a,V 


t  2J    (barns)    9.300  8.733   10.22   9.452   10.60  9.711 
„  -/  (barns)    3.507   3.295   3.856  3.566   3.999   3.664 


cpQ  4  (n/cm2-sec)   4.2764  3.5943   7.8618  6.2336  11.1011  8.4558 
5cp,0.4 

V 

6Dy 

^n 

Mn 

xv 

6V 

a/ 

—  W.  E.  Schilling,  private  communication  (August,  1964). 

—  C.  R.  MacVean,  footnote  2,  p.  141. 


1.190 

1.261 

1.313 

0.989 

0.991 

0.993 

1.269 

1.371 

1.446 

0.953 

0.961 

0.967 

1.208 

1.310 

1.386 

0.980 

0.983 

0.987 

1.241 

1.340 

1.414 
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the  reactor  at  two  different  temperatures  (20°  and  40°C),  the 
temperature  dependent  behavior  of  these  properties  was  pre- 
dicted.  Table  4.4  gives  the  calculated  values  of  the  neutron 
density  disadvantage  factor  at  20°C,  the  logarithmic  slope 
(1/6  )(d&  /dO),  and  the  slope  d&n/dO. 


TABLE  4.4 

TEMPERATURE  DEPENDENCE  OF  6 

n 


W/F  6n  SJ  (i/5n)(d6n/dO)-/        d5n/d0 


in  10"4/°C        in  10"4/°C 


1:1 

1.252 

2:1 

1.351 

3:1 

1.422 

-2.76  -3.46 

-4.46  -6.02 

-5.56  -7.91 


a/ 

—  C.  R.  MacVean,  footnote  2,  pp.  141-143. 


CHAPTER  5 

EXPERIMENTAL  METHODS  AND  RESULTS  FOR  THERMAL  DISADVANTAGE 

FACTORS 
5 . 1   REACTOR  DESCRIPTION 
5.1.1   The  Unit  Cell 

The  lattices  of  the  Cornell  University  Zero  Power 
Reactor  consist  of  an  upper  fuel  rod  support  grid  plate, 
a  lower  fuel  rod  positioning  grid  plate,  structural  support 
members,  the  reactor  tank,  the  fuel  rods,  and  four  control 
rod  units  each  containing  a  cluster  of  three  control  ele- 
ments in  an  equilateral  triangle  of  the  same  pitch  as  the 
rest  of  the  core.   These  control  rod  clusters  may  be 
positioned  in  any  location  of  the  reactor  because  the 
control  elements  have  the  same  outside  diameter  as  a  normal 
fuel  rod.   The  upper  portion  of  each  control  element  con- 
tains the  boron  carbide  poison.   The  lower  portion  is  a  fuel 
follower  with  the  same  dimensions  as  the  normal  fuel  rods. 

The  upper  grid  plate  has  an  insert  which  allows  the 
removal  or  insertion  of  four  fuel  rods  as  a  unit  in  what 
is  called  the  unit  cell.   These  four  fuel  rods  form  a 
parallelogram  of  two  equilateral  triangles  of  the  lattice 
pitch.   For  all  water-to-fuel  ratios,  the  unit  cell  is 
slightly  off-center  with  respect  to  the  upper  grid  plate 
and  the  reactor  tank.   All  cores  studied  were  assembled 

-70- 


Sec.  5.1  71 

with  the  unit  cell  at  the  center  of  the  fuel  loading  to  reduce 
the  radial  flux  gradient  in  the  unit  cell .   With  a  triangular 
(or  hexagonal)  lattice,  one  fuel  rod  must  be  the  exact  center 
of  the  core.   Therefore,  of  the  four  fuel  rods  of  the  unit 
cell,  one  was  the  central  element  of  the  reactor  and  the  other 
three  were  in  the  first  hexagon  out. 

5.1.2   Activation  Rods 

The  original  design  of  the  reactor  facility  included  de- 
cappable  fuel  rods  for  the  purpose  of  introducing  foils  into 
the  fuel  region  of  the  reactor.   These  rods  were  slightly 
longer  than  the  normal  fuel  rods  but  contained  a  normal  48- 
inch  fuel  loading. 

Because  the  full-length  decappable  fuel  rod  is  too  long 
to  fit  into  the  hood  used  for  handling  radioactive  materials, 
and  because  it  is  difficult  to  insert  and  remove  foils  with 
this  length  of  fuel  rod,  a  special  "activation"  rod  was  de- 
signed.  The  eighteen- inch  fuel  section  of  the  activation  rod 
is  supported  by  screw-in  aluminum  rods  which  make  this  special 
rod  the  same  over-all  length  as  a  standard  fuel  element. 
(See  Figure  5.1)   When  detached,  the  fuel  section  is  short 
enough  to  be  handled  easily  in  the  hood.   It  is  easy  to  insert 
and  remove  foils  with  this  rod  because  only  a  few  fuel  pellets 
need  be  handled.   Because  of  the  small  size,  ease  in  handling, 
and  small  number  of  fuel  pellets  handled,  experiments  can  be 
accomplished  with  less  radiation  exposure  to  personnel  than 
was  previously  possible  with  the  longer  rods. 


fl 
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HORIZONTAL    MIDPLANE    OF    FUEL    IN 
THE   REACTOR 


BOTTOM   FUEL   PELLET    IN    THE 
ACTIVATION    ROD 

A   IS    THE  DISTANCE   DOWN   FROM    THE 
REACTOR   HORIZONTAL    FUEL   MIDPLANE 
TO  THE   BOTTOM   FACE   OF   THE    BOTTOM 
FUEL    PELLET  IN   THE  ACTIVATION  RODS. 


W/F 
III 

2!l 
3!l 


0" 
3.5' 
9.0' 


Figure    5.1:      The  Activation    Rod 


Sec.  5.1  -  5.2  73 

The  activation  rods  are  used  for  foil  and  wire  activation 
within  the  fuel  region  of  the  core.   They  can  be  placed  in  any 
normal  fuel  rod  position,  but  for  the  disadvantage  factor  work 
they  were  used  only  in  the  unit  cell.   The  four  fuel  rod 
positions  of  the  unit  cell  were  filled  with  three  activation 
rods  and  one  normal  fuel  rod.   The  standard  loading  was  as 
shown  in  Figure  5.2. 

5.2   FOIL  PREPARATION 
5.2.1   Foil  Types  and  Manufacture 

The  foils  used  for  the  disadvantage  factor  measurement 
were  five  mil- thick  dysprosium,  manganese,  and  vanadium. 
The  dysprosium  was  obtained  from  General  Astrometals  Corpora- 
tion and  consisted  of  5.24%  dysprosium  in  aluminum.   The 
manganese  was  in  the  form  of  an  alloy  called  P-metal  which 
is  about  72%  manganese,  18%  copper,  and  10%  nickel.   The 
P-metal  was  obtained  from  Metals  and  Controls,  Inc.,  and  has 
a  thickness  tolerance  of  ±0.00025  inch.   The  vanadium  was 
obtained  from  H.  Cross  Co. 

The  fuel  region  foils  were  punched  from  stock  with  a 
0 .600±0 .001-inch  diameter  die.   Moderator  region  foils  were 
punched  from  stock  using  a  machined  die  of  appropriate 
shape.   Figure  5.3  shows  the  dimensions  and  shape  of  the 
moderator  foils. 
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A    ACTIVATION    ROD    IN    CENTRAL    FUEL   POSITION 
OF    THE  UNIT    CELL. 

B,C  ACTIVATION    RODS    IN    FIRST     HEXAGON   POSITIONS 
OF   THE   UNIT    CELL. 

D    NORMAL    FUEL  ELEMENT    IN    A    FIRST    HEXAGON 
POSITION    OF    THE  UNIT    CELL. 

E    REMOVABLE    UNIT    CELL    FUEL    ROD     SUPPORT 
PLATE. 

F     NORMAL    FUEL   RODS    IN    NORMAL    LATTICE    POSITION 


Figure  5.2:   Top  View  of  the  Unit  Cell 
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W/F 

A 
(in) 

B 
(in) 

i:  i 

0.675 

0.856 

2: 1 

0.675 

1.030 

3: 1 

0.675 

1.177 

Figure  5.3:   Moderator  Foil  Dimensions 
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5.2.2   Foil  Placement  in  the  Reactor 

To  prevent  fission  product  contamination  of  the  fuel 
foils,  a  special  fuel  pellet  with  an  aluminum  foil  glued 
to  one  end  with  G.E.  Glyptal  was  placed  on  each  side  of  the 
fuel  foils.   See  Figure  5.4.   Ten  fuel  pellets  were  removed 
from  the  activation  rod  and  an  aluminum  lined  fuel  pellet 
was  inserted.   The  manganese  fuel  foil  was  placed  next  to 
the  fuel  pellet  liner.   Pellet  and  foil  were  pushed  to  full 
insertion  with  a  lucite  rod  which  was  marked  for  trans- 
ferring the  foil  position  to  the  outside  of  the  activation 
rod.   Another  lined  pellet  was  inserted,  followed  by  a 
normal  fuel  pellet.   This  loading  sequence  was  continued 
for  dysprosium  and  vanadium  foils  until  each  activation 
rod  had  a  manganese  foil  at  the  center  of  the  rod,  a 
dysprosium  foil  three  inches  below  the  center  of  the  rod, 
and  a  vanadium  foil  another  three  inches  lower.   Figure  5.1 
gives  the  distance  between  the  horizontal  midplane  of  the 
fuel  in  the  reactor  and  the  bottom  edge  of  the  bottom  fuel 
pellet  in  the  activation  rod  for  the  three  water-to-fuel 
ratios  studied. 

Two  moderator  foils  were  secured  to  a  one-eighth-inch 
thick  lucite  moderator  support  plate  with  double-backed 
Scotch  tape.   The  unit  cell  was  assembled  with  three  acti- 
vation rods  and  a  normal  fuel  rod.   Using  the  marked  lucite 
fuel  foil  insertion  rod,  the  lucite  moderator  foil  supports 
were  positioned  at  the  proper  levels  on  the  outside  of  the 
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Figure  5.4:   Fuel  Foil  Placement  in  the  Activation  Rod 
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activation  rods.   The  assembly  was  taped  together  with  Band 
Aid  Air  Vent  adhesive  tape  placed  about  three  inches  above 
the  manganese  foils  and  about  three  inches  lower  than  the 
vanadium  foils.   (This  tape  was  found  to  be  easy  to  use,  to 
keep  its  adhesive  properties  under  water,  and  to  come  cleanly 
off  the  fuel  rods  leaving  no  sticky  residue.)   To  insure 
positive  positioning  of  the  lucite  moderator  support  plates, 
and  to  prevent  loss  of  foils,  waxed  string  was  looped 
around  the  fuel  rods  immediately  above  and  below  the  support 
unit.   Figure  5.5  shows  the  assembly  ready  for  insertion 
into  the  reactor. 

The  assembled  unit  cell  was  then  loaded  into  the  reactor 
and  irradiated  for  a  period  of  about  ten  to  twenty  minutes 
at  flux  levels  corresponding  to  about  one-half  watt  in  power. 
Irradiation  times  and  power  levels  were  adjusted  to  give 
good  count  rates  with  minimum  personnel  exposure  during 
handling  of  the  fuel . 

It  was  necessary  to  remove  the  unit  cell  assembly  as 
rapidly  as  possible  to  get  any  information  at  all  from  the 
short-lived  vanadium  foils  (half  life  of  3.8  minutes). 
After  a  wait  of  about  five  minutes  from  shutdown,  the  acti- 
vation rods  were  removed  to  the  radioactive  handling  hood 
for  disassembly.   The  vanadium  foils  were  removed,  washed 
in  acetone,  and  beta-counted  on  an  automatic  sample  changer. 
About  one  hour  after  shutdown,  the  manganese  and  dysprosium 
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foils  were  removed  from  the  fuel,  washed  in  acetone,  and 
counted  on  the  automatic  sample  changer. 

5.3   COUNTING  TECHNIQUES  --  SOLID  FOILS 

5.3.1  The  Automatic  Counter  System 

38  / 
The  automatic  counting  system  described  by  Rushing — 

39/ 
and  Berg- —  was  used  for  the  solid  foil  counting.   There 

are  24  sample  positions  on  a  round  table  which  pneumatically 
advances  on  electronic  signal.   The  samples  are  held  on 
pedestals  which  are  pushed  up  next  to  a  thin-window  pro- 
portional flow  counter  for  counting.   The  timing  interval, 
table  rotation,  two  scalers,  and  a  printer  are  controlled 
electronically . 

Uniform  counting  geometry  for  the  fuel  foils  was  pro- 
vided by  pedestals  which  were  machined  to  accept  0.600- 
inch  foils  o   The  moderator  foils  were  centered  by  eye  on 
similar  pedestals. 

5.3.2  Relative  Fuel  and  Moderator  Foil  Counter  Efficiencies 

Counter  efficiency  is  slightly  sensitive  to  distance 

40/ 
from  the  center  of  the  pedestal.—   Therefore  an  experi- 


38  / 

—  J.  C.  Rushing,  "Report  of  Work  Done  on  an  Automatic 

Sample  Changer  for  Use  in  the  Cornell  University 

Nuclear  Reactor  Laboratory,"  Department  of  Engineering 

Physics  Senior  Project  (1962). 

39/ 

— '  S.  S.  Berg,  footnote  1,  pp.  101-107. 
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mental  geometric  correction  factor  was  determined  for  this 
effect  as  follows.   Fuel  and  moderator  foils  were  irradiated 

on  a  flux  wheel  in  the  uniform  thermal  flux  of  the  TRIGA 

41/ 
hohlraum—  and  both  shapes  were  counted  on  the  automatic 

counting  system.   After  irradiation  in  a  uniform  flux,  the 

two  shapes  of  foils  should  have  the  same  activation  per  unit 

volume.   The  foils  were  counted  using  the  Geneva  motion 

described  in  Section  5.3.4.   The  ratio  of  fuel  foil  count 

rate  per  unit  volume  to  the  moderator  foil  count  rate  per 

unit  volume  is  the  desired  geometrical  correction  factor. 

These  correction  factors  are  given  in  Table  5.1  for  the  three 

water-to-fuel  ratios  studied. 

TABLE  5.1 

EXPERIMENTALLY  DETERMINED  COUNTER  EFFICIENCY  CORRECTION 

FACTORS 

W/F  1:1  2:1  3:1 

Factor  -1  0.990  1.007  1.014 

a/ 

—  The  statistical  error  is  about  0.004  for  all  factors. 


—  E.  Bo  Fehr,  internal  memo  RM-16,  Cornell  University 
Nuclear  Reactor  Laboratory  (1964)  . 

41/ 

—'    E.  B.  Fehr  and  E.  P.  Whit ted,  internal  memo  RM-15, 

Cornell  University  Reactor  Laboratory  (1964) . 
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5.3.3  Data  Reduction  for  Vanadium  Foils 

The  short  half-life  of  the  vanadium  foils  made  it  necessary 
to  analyze  the  data  differently  than  the  data  from  the  other  two 
foils .   The  vanadium  foils  were  counted  in  a  given  order  on  both 
proportional  counters  of  the  automatic  counting  system.   The 
count  rate  at  the  start  of  counting  was  about  1,000  counts  per 

second,   The  foils  were  counted  down  to  about  20  counts  per 

42/ 
second.   SAT  I,—  a  local  computer  code,  was  used  to  correct 

each  datum  point  back  to  the  activity  on  removal  from  the 

reactor.   These  count  rates  on  removal  were  then  converted  to 

count  rate  per  unit  volume.   The  moderator  foil  count  rates 

were  averaged,  as  were  the  fuel  foil  count  rates.   The  ratio 

of  these  averaged  count  rates  per  unit  volume  for  the  moderator 

and  fuel  foils  was  corrected  by  the  appropriate  geometrical 

correction  factor  of  Table  5.1;  the  result  was  the  desired 

vanadium  disadvantage  factor  5  „ 

5.3.4  Data  Reduction  for  Manganese  and  Dysprosium  Foils 
Manganese  and  dysprosium  have  half-lives  of  2.59  hours 

and  2.42  hours  respectively.   This  is  sufficiently  long  to 
allow  use  of  the  Geneva  motion  counting  technique.   In  this 
method,  the  foils  are  counted  in  a  given  order  once  through. 
The  last  foil  is  counted  twice.   The  remaining  foils  are  then 
counted  in  the  reversed  order.   For  a  total  counting  time 


42/ 

—  W.  E.  Schilling,  internal  memo,  Cornell  University  Reactor 

Laboratory  (1962)  „ 
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sufficiently  short  compared  with  the  half-life  of  the  foils, 
a  simple  arithmetic  time  average  of  two  counts  is  very  close 
to  the  actual  counts  at  the  midpoint  in  time.   For  example, 
for  a  total  counting  time  of  20  minutes,  the  arithmetic 
average  of  the  end  point  counts  is  within  0.3%  of  the  true 
counts  at  the  midpoint  time.   The  total  counting  time  was 
usually  about  10  minutes;  therefore  the  error  is  even  less. 
Statistical  accuracy  of  counting  was  generally  of  the  order 
0.7%  and  thus  the  error  introduced  by  using  the  Geneva 
counting  method  is  much  less  than  the  statistical  error. 

5.4   COUNTING  TECHNIQUE  --  DISSOLVED  FOILS 
The  manganese  foils  were  counted  with  a  "soluble" 
method  in  which  the  foils  were  dissolved  in  acid  and 
aliquots  of  the  solution  were  gamma-counted  in  a  well- 
type  scintillation  counter.   Counting  of  aliquots  of 
solutions  of  foils  provides  two  advantages  over  solid 
foil  counting.   First,  the  moderator  and  fuel  foils  are 
counted  under  identical  counting  geometry.   Thus,  there 
is  no  need  to  determine  a  counter  efficiency  correction 
factor.   Second,  the  process  of  dissolving  the  foil 
averages  the  count  rate  of  the  foil  over  its  volume.   The 
count  rate  of  the  aliquot  is  then  directly  proportional 
to  the  total  count  rate  of  the  foil  itself,  the  constant 
of  proportionality  being  the  fraction  of  the  total  solution 
which  is  being  counted.   For  both  moderator  and  fuel  foils, 
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the  constant  of  proportionality  was  the  same.   The  above  total 
count  rate,  when  divided  by  the  volume  of  the  foil,  is  exactly 
the  spatially  averaged  count  rate  per  unit  volume  needed  to 
determine  the  manganese  disadvantage  factor. 

5.4.1  Foil  Preparation 

After  the  solid  foils  had  been  counted  as  described  in 
Section  5.3.4,  they  were  removed  to  the  fume  hood  in  the 
isotope  handling  room.   Each  foil  was  cut  with  metal  shears 
so  that  it  would  fit  down  the  neck  of  a  ten-milliter 
volumetric  flask  into  which  about  five  milliliters  of 
3N  HNOo  had  previously  been  placed.   Once  the  foil  was  dis- 
solved, the  solution  was  diluted  to  10  milliliters  with  dis- 
tilled water.   Each  solution  was  then  mixed  well. 

Three  milliliters  of  the  solution  were  pipetted  off 
and  placed  into  a  plastic  vial,  and  the  vial  was  capped. 
Each  foil  had  its  own  vial;  each  vial  had  three  milliliters 
of  solution.   Two  pipettes  were  used,  one  for  the  fuel 
toils  and  one  for  the  moderator  foils  to  prevent  mixing  in 
the  pipette  of  solutions  with  different  activities. 

5.4.2  Counting  of  the  Solutions 

The  vials  fit  snugly  into  the  well  of  the  scintillation 
counter,  the  top  surface  of  the  solution  being  below  the  top 
of  the  crystal.   The  geometrical  arrangement  in  the  counter 
was  identical  for  each  solution. 
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Three  counting  methods  were  initially  investigated. 
First,  a  multichannel  analyzer  was  used  to  get  the  gamma 
spectrum  of  the  solutions.   The  area  under  the  peak  was 
determined  by  numerical  integration.   These  total  counts 
were  divided  by  the  foil  weight  to  get  counts  per  unit 
weight,  and  the  ratio  of  the  average  counts  per  unit 
weight  for  moderator  and  fuel  foils  gave  the  experimental 
disadvantage  factor. 

The  second  method  was  to  count  the  solution  on  a 
single  channel  analyzer  in  the  differential  mode,  adjusting 
the  base  line  and  window  width  to  straddle  the  gamma  peak. 
Again,  the  ratio  of  counts  per  unit  weight  for  moderator 
and  fuel  foils  was  the  desired  disadvantage  factor. 

The  third  method  used  an  integral  discriminator  with 
the  base  line  in  the  valley  above  the  Compton  edge.   In 
this  manner,  all  counts  above  the  base  line  were  included 
in  the  count  rate  per  unit  weight. 

All  three  methods  gave  the  same  result  within  experi- 
mental error.   The  third  method  was  chosen  because  of 
simplicity  of  analysis,  greater  availability  of  the 
associated  electronic  equipment,  and  higher  count  rate 
which  resulted  in  greater  statistical  accuracy. 

The  aliquots  were  counted  using  the  Geneva  motion  pre- 
viously described.   Typical  counting  times  were  one  minute 
per  vial;  typical  total  counts  were  of  the  order  50,000. 
The  background  counts  were  subtracted  from  the  total  counts 
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prior  to  averaging  the  counts  per  unit  weight  for  the  fuel 
and  moderator  foils.  A  typical  data  reduction  is  given  in 
Table  5.2. 

5.5   EXPERIMENTAL  RESULTS 
Disadvantage  factors  were  measured  as  previously  de- 
scribed at  various  temperatures  for  three  water-to-fuel 
ratios.   The  experimental  results  are  given  in  Table  5.3 

as  5~  ,  5..  ,  5,,     ,  (manganese  soluble),  and  5T7.   All 
Dy'   Mn'   Mn,sol  v    °  ''  V 

values  except  the  dissolved  foil  experiment  have  been 
corrected  for  geometry  by  the  correction  factors  of  Table 
5.1.   Soluble  manganese  experiments  had  identical  counting 
geometry  and  therefore  needed  no  geometric  correction 
factor . 

Values  for  the  disadvantage  factors  from  Table  5.3 
were  least-squares  fitted  to  a  linear  function  of  tempera- 
ture.  Table  5.4  gives  fitted  experimental  values  and 
theoretical  values  of  the  disadvantage  factors  at  20°C. 
Both  experimental  and  theoretical  values  of  the  slope  of 
the  5  vs  0  curves  are  also  given. 

Temperature  dependent  values  of  8    and  b  ,  are 

plotted  in  Figure  5.6  for  the  three  water-to-fuel  ratios 
investigated.   Figure  5.7  gives  b        vs  temperature  for  the 
3:1-2  core.   Figure  5.8  gives  b      vs  temperature  for  the 
same  core. 
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TABLE  5.2 

DATA  ANALYSIS  FOR  A  TYPICAL  MANGANESE  DISADVANTAGE  FACTOR 
USING  DISSOLVED  FOIL  TECHNIQUE  MEASUREMENT 


Single  Channel  Analyzer  Settings,  Integral  Mode 
Fine  Gain     0.5  Base  Line     17.0  V 

Coarse  Gain   4  Window  Width    0.5V 

Clipping  Time  0.5  (microseconds)    Discriminator   8.5  V 


Foil*7 

Counts/ 

Average 

Average 

Foil 

Counts 

Unit 

Counts/ 

Counts/ 

&Mn,Sol 

Number 

Weight 

per 

Weight/ 

Weight, 

Weight, 

(gm) 

Minute 

Minute 

Moderator 
Foils 

Fuel 
Foils 

119 

0.1749 

21402 

237  719 

310  868 

238  376 

1.304+ 

307 

0.2533 

39897 

308  366 

0.007 

124 

0.1748 

21043 

239  273 

317 

0.2534 

40677 

313  370 

131 

0.1749 

21373 

238  136 

131 

20629 

317 

39083 

124 

20774 

307 

38564 

119 

20527 

a/ 

—  Foils  numbered  in  the  100 's  are  fuel  foils;  those  in  the 

300 's  are  moderator  foils.   These  foils  are  for  the 

3 : 1  core . 

b  / 

—  Counts  per  unit  weight  per  minute  are  the  sum  of  counts 

for  each  foil,  corrected  for  background  of  176  counts 
per  minute,  and  divided  by  the  foil  weight  ° 


Sec.    5=5 


88 


CO 

lo 
W 


LO 


lO 


UO 


uO 


®o° 


C1J 

Q 


fa 


o 
p— I 

o 

o 
m 


O 
+1 
<N 

CO 


O 

-4-1 


oo  r-^ 

r- 

O  O 

o 

o  o 

o 

e          * 

o 

o  o 

o 

+  1    H 

41 

i— 1  1^ 

CO 

LO  v£> 

LP) 

o  o  o  o  o  o 
o  o  o  o  o  o 

o  o  o  o  o  o 

fl       hi       -;      +!      +!       4: 

CO  ^O  vO  m  LO  t-O 


o^ 

a> 

o^ 

o 

o 

o 

o 

o 

o 

ft 

e 

e 

o 

O 

o 

1 

+1 

fl 

..) 

C! 

CM 

m 

I-- 

[-». 

vOOOOO^Ot^ 
^O  m  <t  O  <J"  On 

On  CM  cO  CO  i— I  LO 
CM  CM  CO  <t  vO 


<1-  <t-  <t  <t  <f  <j" 

O^O  vO  vO  O  vO 

>X)  o  co  o>  r^  r-i 

r— I     i — I     i — I 

CO  CO  <f  CO  <f  CO 


CN 

I 


o  o  o 
1—1 1— I  I— I 
o  o  o 

o  o  o 

4-1  4-1  4-1 

LO<t  1— I 

vd  co  r^ 

CM  CM  CM 


r». 

r-» 

o 

o 

o 

o 

0 

© 

o 

o 

fl 

fl 

OO 

co 

r-. 

r-. 

CN 

CM 

CO  CO  O  CO 

CO 

LO  LP)  vO  LT|  v£3 

o  o  o  o 

o 

O  O  O  O  O 

o  o  o  o 

o 

O  O  O  O  O 

o  o  o  o 

o 

o  o  o  o  o 

fl    -rl 

+  1 

1  -fl  +1  +1  +1 

m  m  v£>  o 

o> 

<t  O  CM  i— 1  0> 

00  1^  vO  LO 

v£> 

CO  CO  r— 1   i— 1  O 

CM  CM  CM  CM 

CM 

CO  CO  CO  CO  CO 

CM      UO  LP) 

i-lLOOMNHst- 

mHMOMN  H 
CM  <f  LO  CM  CM 


<f  <r  <r  <t  <t  <r 

ON  CO  L(0  <t  CM  v£) 

i—l   i— I   i— (   r— I   i— I  CM 
LP)  UO  LO   LP)  LO   UO 


CM 


CM 


co  t^-  co  co  r*- 
o  o  o  o  o 
o  o  o  o  o 

o  o  o  o  o 

+1    +1     fl     f!     fl 

<t  LO  0>  CM  O 
CO  CO  CM  CO  CO 


o  o  o  o  o 
o  o  o  o  o 

o  o  o  o  o 

fl  fl  +1  +1  fl 

<f  CM  CM  LO  CO 

CO  CO  O  O  i—l 
CO  CO  CO  CO  CO 


o 
1—1 

o 

lo  in  vD  o  rv 
o  o  o  o  o 
o  o  o  o  o 

o 

o  o  o  o  o 

f! 
00 
00 
CM 

O  <j"  rH  CO  CO 

oo  r^.  <t  <j-  lo 

CO  CO  CO  CO  CO 

CM  CM 

CM^OCOCO 

<)■  o  o  r-^  o-n 
CM  co  co  <)- 


^O  ^£>  vQ  *X>  v£> 

co  r-»  oo  r^  vd 


r^  r^.  r-  r^  r^- 


CM 


ro 


Sec.    5.5 


89 


> 


UO 

w 


o 

Em 

c 

H 

fa 

Q 
O 

s 

fa 

o 

H 

<d 

p^ 

o 

Z 
H 

en 

CO 


o 
fa 


a 
X 

a; 


O 
H 

o 

H 

o 

< 
fa 

o 

H 

H 
O 
fa 


O 

o 


--I      C 
CO 


oj 


Q 


> 


Q 


o 
m 


-± 

■—I 
o 

CM 


cm 


m 


o 

CM 

CO 


xD 

CO 


LO 

<t 

CO 

o 

i— I 

a> 
o 


o 
t— I 

uO 
O 


O^ 
CA 

o 

Cxi 

LO 


xD 

CO 

r-i 

ur\ 
O 


oo 

LO 

O 


uO 

I— I 

^D 
<t 
O 


o 

o 

CO 

o 


OO 
CO 

o 


O 


ci 
r^ 
CN 
O 


CO 

o 


vO 
uO 


LO 

CM 

CO 
uO 


uO 
O 


CI 

O 


Ox 
00 
00 
O 


CA 

oo 

UO 
O 


00 

o 


o^ 


uO 


LO 

o^ 

o^ 


xD 
OA 


o 
o^ 


<t 
o^ 


CM 

o> 


o 

CO 
CA 


CO 
CA 

o^ 


% 


•r-l 

•H 

•H 

i 

o 

o 

O 

P4 

V 

> 

> 

TJ 

■~^ 

6 

M-J 

CO 

^H 

> 

> 

X3 

r-l   i— 1 

•H 

cti 

I    o 

O 

ii 

pH 

uO 

M-l  U-i 

0 

+ 

> 

p 

i 

to 

r-l 

pd 


-90- 


o 

o 


<D 


^j 


m 


'D 
U 
D 
4J 

03 

:-. 

a. 

E 

QJ 

H 

W 
> 

5-i 

C 
U 

u 

Pi-! 

0) 

SO 

n3 
U 

c 

> 

r. 
Q 
(U 

X 

y 

C 

bC 

c 
ta 

2C 


m 


CO 

3 


o 


00 

(0 

<*- 

C\J 

O 

<M 

C\J 

C\J 

OJ 

OJ 

to 


UW 


-91- 


CO 


1.40 

1.38 

1.36 
o 

1.34 

1.32 

1.30 

1.28 

W/F  =  3:i-2 


10 


20 


30 


e  (in°o 


40 


50 


Figure  5.7:   Dysprosium  Disadvantage  Factors  vs  Temperature;  3:1-2 
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Figure  5.8:   Vanadium  Disadvantage  Factors  vs  Temperature;  3:1-2 
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5,6   DISCUSSION  OF  RESULTS 
5.6.1   Foil  Construction  and  Positioning  Errors 

Moderator  foils  were  punched  from  stock  using  a  machined 
die  of  appropriate  shape.   The  radius  of  the  circular  cuts 
taken  in  the  corners  of  the  triangles  was  that  of  a  standard 
drill  size  (0.6875  inch)  rather  than  the  exact  radius  of  the 
outside  of  the  fuel  rods  (0.666  inch).   This  gave  all  the 
moderator  foils  a  slightly  smaller  volume  than  the  actual 
moderator  volume  of  interest. 

An  estimate  of  the  error  associated  with  inexactly  cut 
foils  can  be  made  as  follows.   First  it  is  assumed  that  the 
foil  thickness  is  uniform,   A  second  assumption  is  that  the 
only  error  in  foil  construction  will  be  in  making  the  modera- 
tor foil.   This  is  justified  somewhat  in  that  the  fuel  foil 
is  made  from  a  die  with  less  than  0 ,  570  tolerance.   The 
following  equation  is  written  for  the  correct  experimental 
disadvantage  factor  for  a  foil  with  1/v  absorption  cross 
section: 

„  /T?m        T7m  .    ■r=m  TTm   N  /T7m 

77m  (N^    .-.VV    . -,    +  N   .pV    p)/V 

<-  C  foil    foil  nf  nfy  ,r    -,  >. 

correct  ci  ,^-r        yr        ^ /yr 

UNfoilvfoil;/v 

where  the  subscript  "foil"  indicates  a  property  of  the  foil, 
and  the  subscript  nf  indicates  a  property  of  the  missing 
portion  of  the  foil,  or  the  "non-foil."   All  other  symbols 
have  been  previously  defined,   The  second  assumption  above 
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f    f 
sets  V   =  V r    .  -,  .   Equation  (5  1)  can  then  be  rewritten  as 

5        =6      [(V?  .J\Jm)    +  flf^vVfi?  -iVm")l      (5.2) 
correct    meas   v  foil       v  nf  nf   foil  _j      v    y 

where  5     is  the  previously  reported  disadvantage  factor 

ruG  Si  S 

as  given  in  Table  5  3   Assuming  that  the  missing  section 
of  the  foil  is  from  the  circular  region  near  the  fuel 

cladding,  the  neutron  density  cladding  disadvantage  factors 

a  3  /  ="™m     /~Hm 

of  MacVean—      can  be   used    to    approximate   N    r/Nr    .  n  .      This 

rr  nf   foil 

will  result  in  an  overestimate  of  the  correction  factor. 
Let 

R  =  vm/Vfoii  (5-3) 

Then,  Eq,  (5.2)  becomes 

o      _         R"1[l+(R-l)/&  ]  (5.4) 

correct  —  meas  7   c J  v    y 

where  5   is  the  cladding  disadvantage  factor 

r 

Typical  values  of  V c    . , /Vr  .  i  were  taken  from  foils 
/r  foil   foil 

used  in  the  experiments.   Table  5  5  gives  the  values  of 
the  corrections  for  the  three  foils  and  three  water-to- 
fuel  ratios  used    It  can  be  seen  that  the  correction  is 
small,  and  generally  within  the  accuracy  of  the  experiment-, 
Because  it  is  a  small  correction,  and  because  it  is  not  an 
accurately  determined  correction,  it  has  not  been  applied 
to  the  experimentally  measured  disadvantage  factors 


-=-'    C.    Ro  MacVean,  footnote  2,  p,  141 
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Errors  resulting  from  improper  vertical  positioning  of 
the  moderator  foils  in  the  experimental  assembly  are 
expected  to  be  negligible,   The  flux  change  in  the  axial 
direction  at  the  foil  position  in  the  unit  cell  is  of  the 
order  of  one  or  two  per  cent  per  inch.   Positioning  accuracy 
is  estimated  to  be  about  0,1  inch.   The  resulting  error  is 
then  seen  to  be  about  0,2?o  in  count  rate, 

A  small  error  may  be  introduced  by  the  presence  of 
the  lucite  support  plate  for  the  moderator  foil.   This 
plate  is  one-eighth-inch  thick  and  of  approximately  the 
same  shape  as  the  moderator  foils.   It  is  suggested  that 
further  experimental  work  be  done  with  different  thicknesses 
of  lucite  to  investigate  the  magnitude  of  this  effect, 

The  error  associated  with  the  finite  thickness  of  the 
foils  has  not  been  investigated  experimentally.   It  was 
assumed  in  Chapter  4  that  the  self -shielding  factors  for 
the  fuel  and  moderator  foils  were  equal,   In  conjunction 
with  the  investigation  of  lucite  thickness  effect,  it  is 
suggested  that  the  disadvantage  factor  be  measured  as  a 
function  of  detector  foil  thickness. 

The  process  of  dissolving  the  foils  for  the  "soluble" 
disadvantage  factor  measurement  is  not  expected  to  con- 
tribute any  error  to  the  measurement,   The  chemical  handling 
technique  employed  standard  volumetric  flasks  and  volumetric 
pipettes.   The  pipettes  are  estimated  to  be  accurate  to 
about  0.25%  and  the  flasks  to  about  0 . 1% -   The  error  should 
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be  uniform  for  each  aliquot,  and  therefore  cancel  from  the 
ratio  which  is  the  disadvantage  factor. 

Weighing  of  the  foils  was  accurate  to  0.1  mg  which 
represents  a  maximum  error  of  less  than  0.3%  for  moderator 
foils  of  the  1:1  core.   For  all  other  foils  the  error  is 
less  than  0.2%. 

5.6.2   Foil  Induced  Flux  Perturbations 

There  are  two  perturbations  associated  with  inserting 
foils  into  the  reactor  for  measurement.   One  effect  is  the 
depression  of  the  flux  by  the  foil  itself  as  in  the  case 
of  an  absorber  being  placed  into  the  low  absorption  region 
of  the  moderator.   This  would  tend  to  reduce  the  measured 
value  of  the  flux  in  the  moderator.   This  effect  could  be 
determined  by  measuring  the  foil  activity  per  unit  weight 
as  a  function  of  foil  thickness  and  extrapolating  to  zero 
foil  thickness. 

The  second  effect  is  present  in  the  fuel  foil.   Here, 
a  highly  absorbing  region  is  replaced  by  a  lightly  absorbing 
foil.   Neutrons  may  stream  into  the  foil  and  the  average 
flux  in  the  foil  would  then  be  higher  than  that  in  a 
corresponding  fuel  region.   This  effect  could  also  be  deter- 
mined by  extrapolation  to  zero  thickness  foils. 

Both  of  these  effects  tend  to  make  the  measured  dis- 
advantage factor  lower  than  the  ratio  of  the  neutron  den- 
sities appropriate  to  the  two  regions  under  investigation. 
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An  estimate  of  the  first  correction,  called  the 
outer  flux  perturbation,  was  made  using  the  procedure 
outlined  in  ANL  5800.- —   A  correction  factor,  defined 
as  the  ratio  of  the  activation  in  the  perturbed  flux 
to  activation  in  the  unperturbed  flux  was  calculated 
for  the  three  foil  materials  and  for  water- to- fuel 
ratios  of  1:1  and  3:1  using  the  approximations  given 
to  the  transport  theory  calculations  of  Ritchie  and 
Eldridge.   This  correction  is  for  a  coin  shaped  foil, 
but  the  moderator  foils  used  here  were  of  the  shape 
shown  in  Figure  5.3.   An  equivalent  radius  for  these 
foils  was  taken  to  be  0.5  cm  and  1.0  cm  for  the  1:1 
and  3:1  cores  respectively.   The  correction  was  then 
estimated  to  be:  0.3-0.8%  for  dysprosium,  0.6-1. 6% 
for  manganese,  and  0.1-0.47o  for  vanadium  where  the 
lower  values  are  for  the  1:1  core.   The  correction 
increases  the  measured  disadvantage  factors  by  this 
amount . 

Because  of  the  smallness  of  the  corrections,  and 
the  uncertainty  in  their  accuracy,  they  have  not  been 
applied  to  the  measured  disadvantage  factors.   It  is 
recommended,  rather,  that  an  attempt  be  made  to  measure 
them  experimentally. 
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5.6.3   Flux  Perturbations  Due  to  the  Activation  Rods 

Axial  flux  scans  taken  on  the  1:1  core  showed  a  large 
flux  peaking  near  the  ends  of  the  fueled  section  of  the 
activation  rods.   This  peak  is  a  result  of  replacing  the 
fuel  region  of  a  normal  fuel  rod  with  the  aluminum  support 
pieces  of  the  activation  rod.   Peaking  was  most  pronounced 
in  the  1:1  core  because  of  undermoderation  for  this  water- 
to-  fuel  ratio.   The  peaking  may  have  been  sufficient  to 
affect  the  vanadium  and  possibly  the  dysprosium  disadvantage 
factors  in  the  1:1  core.   The  manganese  disadvantage  factors 
are  not  expected  to  be  influenced  by  this  peaking  because 
these  foils  were  in  the  center  of  the  fueled  section  of  the 
activation  rod. 

Peaking  was  less  pronounced  in  the  2:1  and  3:1  cores. 
In  addition,  the  fueled  section  of  the  activation  rods  was 
lowered  in  successive  runs  as  indicated  in  Figure  5.1. 
The  first  change  was  made  so  that  the  vanadium  foils  were 
in  the  center  of  the  core.   This  was  changed  again  to  the 
"symmetric"  activation  rod  wherein  the  center  of  the  fueled 
section  of  the  activation  rod  is  at  the  center  of  the  fuel 
region  of  the  core. 

The  effect  of  the  peaking  may  result  in  some  inaccuracy 
in  the  measured  value  of  the  disadvantage  factors  for 
vanadium  and  possibly  dysprosium,  but  the  temperature  de- 
pendence of  the  disadvantage  factors  should  not  be  affected. 
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The  only  way  to  remove  this  effect  is  to  resume  use  of  the 
full  length  decappable  fuel  rods .   A  minimization  of  the 
effect  can  be  achieved  by  using  only  one  set  of  foils  and 
placing  them  at  the  center  of  the  fueled  region  of  the 
activation  rod. 

5.6.4   Comparison  of  Experiment  and  Theory 

Table  5.5  compares  theory  and  experiment  and  shows 
that  (1)  dysprosium  disadvantage  factors  lie  between  4 
and  8%  lower  than  theory,  (2)  manganese  disadvantage 
factors  lie  between  2.4  and  4.570  lower  than  theory,  and 
(3)  vanadium  disadvantage  factors  lie  about  4%  lower  than 
theory.   However,  all  measured  values  of  the  temperature 
dependence  of  the  disadvantage  factors  agree  with  the 
theoretical  values  within  the  limits  of  experimental  error. 

5*7   CONCLUSIONS 
The  disadvantage  factor  was  measured  for  three 
water-to-fuel  ratios  using  foils  of  three  different 
materials.   Values  of  the  measured  disadvantage  factors 
are  consistently  a  few  per  cent  lower  than  the  calculated 
values.   Further  refinements  in  the  experimental  technique, 
specifically  utilization  of  different  thickness  foils  and 
moderator  support  plates,  could  reasonably  be  expected  to 
bring  the  experimental  values  of  the  disadvantage  factor 
into  closer  agreement  with  the  calculated  values. 
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The  temperature  dependence  of  the  disadvantage  factor 
was  measured  for  three  water-to-fuel  ratios.   Measured 
values  of  the  temperature  dependence  agree  within  experi- 
mental errors  with  the  theoretical  values  as  calculated  by 
MacVean . — J 
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CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 
6.1   ISOTHERMAL  TEMPERATURE  COEFFICIENT  OF  REACTIVITY 
6.1.1   Summary 

The  isothermal  temperature  coefficient  a  is  defined 
as  dp/d0  where  p  is  reactivity  and  0  is  the  temperature 
of  all  parts  of  the  core.   The  isothermal  temperature  co- 
efficient was  measured  by  differencing,  at  a  series  of 

closely  spaced  values  of  0,  values  of  p  computed  from 

45/ 
observed  clean-core  periods  using  a  local  code  RHOFIT >— 

Temperatures  were  determined  by  averaging  readings  from 

four  thermocouples  distributed  through  the  moderator; 

auxiliary  tests  showed  the  temperature  could  be  made 

uniform  within  iOol°C. 

A  preliminary  calculation  of  a  was  carried  out 

using  two  group,  two  region,  first  order  perturbation 

46/ 
theory  as  outlined  by  Larrimore .— —   A  computer  program 

calculated  the  fluxes  and  adjoint  fluxes  and  used  them 

in  the  standard  perturbation  theory  reactivity  equations 

to  find  the  isothermal  temperature  coefficient. 

Dependence  of  a  on  temperature  can  be  ascribed  to  two 

contributions:  a  rapidly  varying  one  proportional  to 

the  volumetric  coefficient  of  thermal  expansion  of  water, 
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—  W.  E.  Schilling  and  D,  H.  Bryce,  private  communication 

(March,  1964). 
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and  a  second  one,  nearly  constant  over  the  investigated 
temperature  range,  which  is  due  to  temperature  dependence 
of  microscopic  cross  sections  and  other  small  temperature 
dependent  effects.   These  considerations  lead  to  an  equation 
for  a  of  the  form 

a  =  A+BpR  Q  ,  (6.1) 

where  p„  ~  is  the  volumetric  coefficient  of  thermal  expan- 

H2° 
sion  for  water. 

Experimental  temperature  coefficients  were  determined 

for  nominal  water-to-fuel  ratios  of  1:1,  2:1,  3:1,  and 

4:1.   Statistical  errors  in  values  of  a  are  generally  less 

than  107o.   For  the  cores  studied,  and  over  the  temperature 

ranges  used  (1:1,  5-70°C;  2:1,  5-78°C;  3:1,  5-62°C; 

4:1,  7-45°C),  the  isothermal  temperature  coefficient  is 

well  represented  by  Eq .  (6.1)  where  A  and  B  are  different 

constants  for  each  core.   A  comparison  of  theoretical 

values  of  A  and  B  with  experimental  values  found  from  a 

least-squares  fit  to  Eq .  (6.1)  shows  a  need  for  more 

sophisticated  calculations  to  obtain  better  agreement  with 

experiment . 

6.1.2   Conclusions 

Temperature  dependent  excess  reactivities  for  various 
fuel  loadings  were  determined  for  four  water- to- fuel  ratios 
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of  the  Cornell  University  Zero  Power  Reactor  with  sufficient 
accuracy  to  allow  planning  of  further  temperature  dependence 
studies . 

The  isothermal  temperature  coefficient  of  reactivity 
was  measured  for  four  water- to- fuel  ratios  with  an  accuracy 
which  is  comparable  to  that  previously  reported  in  the 
literature.   The  predicted  and  verified  linear  dependence 
of  a  on  ~„  n  had  not  been  reported  previously. 

Since  the  experimental  results  are  based  on  period 
measurements  on  a  clean  core,  and  since  the  calculations 
were  done  with  a  simple  two  group  model,  the  disagreement 
between  calculated  and  measured  values  of  A  and  B  indicates 
the  need  for  a  refinement  in  the  calculations . 

6.2   TEMPERATURE  DEPENDENCE  OF  THERMAL  DISADVANTAGE  FACTORS 
6.2.1   Summary 

The  neutron  density  disadvantage  factor  5   and  the 
experimentally  measured  disadvantage  factor  5     are 

r  J  &  exp 

defined  by  the  ratios  5   =  l^/Nf   and  5    =  (fVc^  where  N 
J  n  exp 

and  C  are  respectively  the  spatially  averaged  neutron 
density  and  foil  count  rate  per  unit  volume;  superscripts 
m  and  f  indicate  respectively  the  moderator  and  fuel 

regions  of  the  reactor. 

i  47/ 

MacVean  s  calculation  of  thermal  lattice  parameters — 

gives  the  temperature  dependence  of  the  disadvantage  factor. 
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Since  the  calculated  changes  are  small,  highly  accurate 
experiments  were  required  for  comparison.   Experimental 
disadvantage  factors  were  measured  for  three  water-to- 
fuel  ratios  over  the  temperature  range  5-65°C. 

An  integral  foil  counting  method  was  used  employing 
a  0 .600±0 .001-inch  diameter  foil  in  the  fuel  and  an 
appropriately  shaped  moderator  foil.   Five-mil  thick 
foils  of  dysprosium  (5.247Q  dysprosium  in  aluminum), 
manganese  (P-metal),  and  vanadium  were  chosen  because 
their  cross  sections  are  nearly  1/v.   Two  counting 
methods  were  used.   One  was  to  dissolve  the  manganese  foils 
in  nitric  acid  and  gamma-count  aliquots  of  the  solutions 
in  a  well- type  scintillation  counter.   The  other  method 
was  to  beta-count  solid  foils  in  a  proportional  flow 
counter.   Results  of  the  two  methods  agree  within 
statistical  errors. 

Experimental  disadvantage  factors  were  least-squares 
fitted  to  linear  functions  of  temperature.   Agreement 
between  theory  and  experiment  for  the  temperature  depend- 
ence of  the  disadvantage  factors  is  within  the  experimental 
error  for  all  cores  and  foils,  but  the  disadvantage  factors 
themselves  are  a  few  per  cent  lower  than  the  calculated 

disadvantage  factors  --  a  result  also  found  in  the  Brook- 

48  / 
haven  and  Westinghouse  work  reported  by  Honeck. — 
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6.2.2   Conclusions 

Experimental  determination  of  the  thermal  disadvantage 
factor  was  accomplished  quite  accurately  using  the  technique 
of  dissolving  the  foils  and  gamma -counting  aliquots  of  the 
solutions.   Beta-counting  of  the  solid  foils  gave  satis- 
factory results  providing  that  measured  values  were  corrected 
for  the  different  geometrical  counting  efficiencies  of 
different  shapes  of  foils.   The  error  for  both  methods  was 
less  than  one  per  cent,  an  improvement  on  previously 
reported  measurements . 

Theory  and  experiment  disagree  by  a  few  per  cent  on 
the  value  of  the  disadvantage  factor.   This  has  been  a 
common  result,  as  mentioned  above.   Several  factors  may  con- 
tribute to  this  disagreement:   (1)  flux  perturbations  due 
to  finite  thickness  foils,  (2)  flux  perturbations  due  to 
the  lucite  moderator  foil  support  plate,  (3)  flux  pertur- 
bations due  to  the  use  of  activation  rods  with  a  shorter 
fuel  region  than  normal  fuel  rods,  (4)  inexactly  shaped 
moderator  foils,  and  (5)  a  small  error  in  the  calculation 
of  the  disadvantage  factor.   The  first  two  possibilities 
could  be  studied  by  using  different  thicknesses  of  foils 
and  of  moderator  foil  support  plates.   The  third  could  be 
eliminated  by  using  full  length  decappable  rods,  or 
minimized  by  placing  foils  only  at  the  center  of  the  fueled 
section  of  the  activation  rod.   A  more  expensive  die  could 
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be  made  for  punching  the  moderator  foils  to  eliminate  the 
fourth  error  source.   Further  study  could  be  made  on  the 
calculational  techniques  for  the  disadvantage  factor; 
however  the  estimated  accuracy  of  27a   for  the  calculated 

values  of  the  neutron  density  disadvantage  factors  by 

49/ 
MacVean —  indicate  that  the  calculations  probably  would 

not  change  much  with  slight  refinements  to  the  code. 

Agreement  between  theory  and  experiment  for  the 

temperature  dependence  of  the  disadvantage  factor  is 

within  the  experimental  error;  however  the  experimental 

error  in  these  results  varies  from  about  207o  to  about 

907o   because  small  differences  were  taken  between  relatively 

large  numbers  to  obtain  the  temperature  dependence. 

6.3   PROPOSALS  FOR  FURTHER  STUDY 
When  the  perturbation  theory  calculations  of  the 
temperature  coefficient  are  brought  into  agreement  with 
the  experiments,  these  same  calculational  techniques 
can  be  used  to  predict  the  effect  of  adding  soluble 
poisons  to  the  moderator.   It  would  then  be  of  interest 
to  measure  the  temperature  coefficient  of  reactivity 
as  a  function  of  the  distributed  poisons. 

As  indicated  earlier,  it  would  be  of  great  interest 
to  investigate  the  effect  that  changing  the  foil  thick- 
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nesses  would  have  on  the  disadvantage  factor  measurements. 
It  might  be  of  further  interest  to  investigate  the  dis- 
advantage factor,  and  possibly  its  temperature  variation, 
using  resonance  foils  to  study  different  energy  regions. 

In  conjunction  with  the  above  proposals,  it  would 
also  be  interesting  to  vary  the  thickness  of  the  lucite 
moderator  foil  support  plate  and  determine  its  effect  on 
the  measured  value  of  the  disadvantage  factor.   Another 
study  could  be  made  using  polyethylene  support  plates. 

Further  into  the  future,  it  might  be  of  interest  to 
measure  these  same  parameters  on  cores  which  have  different 
geometrical  shapes.   For  example  a  square  or  long,  thin 
lattice  could  be  set  up  instead  of  a  cylinder.   The 
ability  to  properly  predict  the  temperature  coefficient 
for  these  geometries  would  be  a  measure  of  the  accuracy 
of  calculating  the  reflector  contribution  to  the  isothermal 
temperature  coefficient  of  reactivity . 


APPENDIX  A 

LATTICE  PARAMETERS  FOR  THE  CORNELL  UNIVERSITY  ZERO  POWER 

REACTOR 

TABLE  A.l 
CORNELL  UNIVERSITY  ZERO  POWER  REACTOR  FUEL  ELEMENT 

PARAMETERS  -' 


Fuel  Material 

Fuel  Rod  Diameter,  in: 

Fuel  Density,  gm/cc: 


Fuel  Enrichment,  Weight  Percent  U 

235 
Fuel  Enrichment,  Atom  Percent  U 

Oxygen  to  Uranium  Number  Ratio: 

Fuel  Molecular  Weight: 

Cladding  Material: 

Cladding  Density,  gm/cc: 

Cladding  Inner  Diameter,  in: 

Cladding  Thickness,  in: 

Fuel  Element  Diameter,  in: 


235 


Uranium  Dioxide 

0.600±0.003 

10.381+0.52% 

2.070±0.10 

2,096+0.90 

2.01:1 

270.26(physical) 

6061-76  aluminum 

2.70 

0.610+0.002 

0.028+0.003 

0.666+.006 
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APPENDIX  B 
ISOTHERMAL  TEMPERATURE  COEFFICIENT  DATA  REDUCTION 

Tables  B»l  -  B.4  give  the  following  information  in 

tabular  form:,  date  of  the  experiment,  reference  run 

number,  number  of  fuel  rods  N,  temperature  0  in  °C, 

period  T  in  seconds,  excess  reactivity  p       ,  reference 
v  '  y    "excess' 

numbers  of  the  two  runs  used  to  calculate  Ap/A0,  Ap ,  A0, 
Ap/A0  =  a,  the  average  temperature  for  the  two  runs  0, 
and  the  volumetric  coefficient  of  thermal  expansion  of 
water  at  the  average  temperature  |3U  r» .   Figure  B.l  is 
a  curve  of  (3„  ~  vs  0  obtained  by  differencing  values  of 
relative  volume  given  in  the  Handbook  of  Physics  and 
Chemistry .—— 
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